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Integrodifferential- und Integralgleichungen für das Gleichgewicht 
dünner elastischer Schalen”) 


f Von N. A. KILTSCHEWSKIJ 3 


Im Folgenden wird eine auf der Anwendung des Reziprozitätssatzes beruhende Methode zur Aufstellung 

von Integrodifferential- und Integralgleichungen für das Gleichgewicht dünner elastischer Schalen in den 

- Verschiebungen dargelegt. Sie ist eine Modifikation der bekannten Methode von Somigliana [5] und unter- 

scheidet sich von dieser durch die spezielle Wahl eines Systems von „‚Hilfsverschiebungen‘‘, die es erlauben, 

die Verrückung eines Punktes der Schalenmittelfläche als Summe aus der Verschiebung eines entsprechenden 

Punktes einer mittleren Plattenebene — als eines ‚‚Bildes‘‘ der mittleren Schalenfläche — und einer zu- 

sätzlichen Verschiebung auszudrücken, die vor allem von der Krümmung der Schalenmittelfläche abhängt. 

Diese Methode, die vom Verfasser bereits vor zehn Jahren auf zylindrische Schalen angewandt wurde [4c], 
wird jetzt auf Schalen mit beliebig gestalteter Mittelfläche ausgedehnt. 


a a nn lin 
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The paper describes a method, based on the reciprocity theorem, for obtaining the integro-differential and 
integral equations for the displacements pertaining to the state of equilibrium of thin elastic shells. This 
method is a modified version of the well-known method due to Somigliana [5] from which it differs by the 
special choice of a system of “auzxiliary displacements”. Due to this choice, the displacement of any point on 
& the medium surface of the shell may be expressed as the sum of the displacement of the corresponding point 
r on the medium surface of a plate — called the “image surface’ of the medium shell surface — and an addi- 
tional displacement depending mainly on the curvature of the medium surface of the shell. T'his method, used 
by the author 10 years ago for cylindrical shells [40], is now extended to apply to shells with a medium surface - 
of any shape whatsoever. 
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Hoske MU31ATAeTCA MEeTON COCTABAIECHHA HHTETPONHUPPEPEeHLMAAIBHEIX H UHTETPAJIBHEIX YPaB- r 
HeHNH PABHOBECHA TOHKUX YIPYTUX 000J104YeR B IIEPeMeleHNHAX, OCHOBAHHBIH HA IPHMeHeHNH e 
TeopeMbI O B3AHMHOCTH PAa60T. ITOT METON ABJIAHETCH BHNONSMEHEHHEM U3BECTHOTO MeToNa 
CoMuAbAHa [5], OTIMYaAch OT MeToNa CoMHNBAHA CIHEIMANBHBIM BbIOOPOM CHCTEMBI ‚‚BCIIOMOTA- 
TeJIbHBIX‘‘ IepeMelleHHi, II03BOJIAIOIUHM IPENCTABHTb IlepeMelleHne TOYEH CPeAMHHOÜ IIO- 
BEPXHOCTU O60JIOYKH CYMMOÜ IIepeMeINeHHA COOTBETCTByIINeH TOYKU CPeMHHHOHÄ ILIOCKOCTH 
IIJACTUHRU, — ‚‚KApTbl‘‘“ CPeNHHHOÜ IIOBEPXHOCTH O00JI0YKH, U MONONHHTEJIBHOTO TTepeMele- 
HUA, 3AaBHUCHIEeTO, B YACTHOCTH, OT KPHBU3HbI CPelHMHHONH HOBEPXHOCTH O00010YRU. HR IMm- 
AIMHNPUYeCKUM 060JIOYKAaM 3TOT METOAX ÖbILJI IIPHMEeHEH ABTOPOM CBblllle NecATu IeT Hasan [4c]. 


B HacToameü pa6oTe yKasabIBaAeTcH PpacııpocTpaHeHHe ITOTO MEeTona HA O6010YKH C IIPON3- 
BOAbBHONU PoPpMOH CpenHuHHOoH IOBEPXHOCTH. 


$ 1. Einleitung. Integralgleiehung für das Gleichgewicht eines Kreisbogens 


Wir zeigen zunächst an dem einfachsten Beispiel einer eindimensionalen Aufgabe gewisse 
Besonderheiten der vorzuschlagenden Methode, um sie dann auf zweidimensionale Probleme des 
Gleichgewichts elastischer Schalen zu erweitern. 

Wir betrachten das Gleichgewicht eines Kreisbogens, auf den eine zu seiner Achse senkrecht 


gerichtete Einzelkraft P vom Betrage 1 wirkt (Bild 1). Er soll als dünner, an beiden Enden 
-befestigter Stab behandelt werden. Die Länge der nicht-deformierten Stabachse sei I. Als 
Lagekoordinate für die Punkte der Stabachse benutzen wir die von A aus zählende Bogenlänge s. 
Für den Angriffspunkt der Kraft P habe sie den Wert 7. Den Punkten des Kreisbogens lassen 
wir jetzt die Punkte der Achse A,B, eines Balkens der Länge I! als Bildpunkte der Punkte der 
Bogenachse entsprechen (Bild 2). Die Lage der Punkte auf der nicht-deformierten Balkenachse 
beschreiben wir durch die von A, aus zählende Längenkoordinate s. Die Koordinatenmaßstäbe 
auf AB und A,B, seien dieselben. 

Der Balken sei an beiden Enden gestützt, und es greife an ihm an der Stelle £ eine zur 


nicht-deformierten Balkenachse A,B, senkrechte Einzellast P vom Betrage 1 an. Mittels des 
Reziprozitätssatzes stellen wir die Beziehungen zwischen den radialen Verschiebungen u(s, 7) 
der Bogenachse und den Durchbiegungen y(s, &) der Balkenachse auf. Die Biegesteifigkeit E J längs 


*) Nach einem am 19. 12. 57 im Institut für Mechanik der Akademie der Wissenschaften der UdSSR 


gehaltenen Vortrag. 
Kl 
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des Bogens werde dabei als konstant und gleich derjenigen des Balkens vorausgesetzt. Die Ver- 
rückungen in Richtung der Bogenachse werden vernachlässigt. 


Die elementare Biegetheorie des Balkens liefert: 


v9 = 9 erı-e— 9) DE rue Be: 


y(s, & a @sI—&—8) Gänge + (1b). 


Bild 2 
Die Drehwinkel des Balkens an den Enden bestimmen sich zu 
AA, FÜ HEUTE 
o(0, I) ee 6 E I ’ oll, &) Den 6 m EA e - . . (1.2). 


Wir wollen jetzt y(s, &) als radiale Verschiebung des Punktes der Bogenachse be- \ 
trachten. Dementsprechend sollen sich auch die Drehwinkel (0, &) und O(l, £) an den Bogen- I 
enden A und B aus (1.2) ergeben. Am Bogen wird die Verschiebung y(s, &) durch die Einzellast N 
und eine zusätzliche, von der Krümmung abhängige kontinuierlich verteilte Radiallast hervor- ; 
gerufen. Diese Zusatzbelastung g(s, &) läßt sich nach der Theorie der dünnen Stäbe aus den Gleich- ö 
gewichtsbeziehungen [5] i 
dy 1dy 
Elta) 0 sonen.) N 
t 
bestimmen. Auf Grund von (1.la, b) finden wir: 
t 
1 = 
ee j 
a?l 

sd) = TEN. j 

al S (l s) (s >&) 


Offensichtlich ist g(s, &) = 4(&; $). 


Wenden wir den Reziprozitätssatz an, so erhalten wir als erstes auf den Bogen wirkendes 
Belastungssystem die in C(n) angreifende Einzelkraft P (vom Betrage 1) und die durch sie hervor- 
gerufenen Reaktionskräfte R(0,n), R(l,n) und Reaktionsmomente M(0,n), M(l,n) (Bild 1). 
Die dadurch hervorgerufenen radialen Verschiebungen seien u(s,n). Das zweite Belastungs- 
system, das Hilfssystem, besteht aus einer Einheitskraft, die in dem dem Punkte D(£) des Balkens 
entsprechenden Punkt D(£) des Bogens angreift, aus der kontinuierlichen radialen Belastung 
g(s, 8), die sich aus (1.3) bestimmt, den Reaktionskräften 0(0, &), (Ol, E) und den Reaktions- 
momenten L(0, &), L(l, &). Die durch dieses System verursachten radialen Verrückungen sind 


y(s, &). 


Bei Anwendung des Reziprozitätssatzes auf diese Belastungssysteme finden wir 


u(&, 7) = y 9 — f ds 9 ucs, 7) ds + M(0, 7) 00,9) + Mn) OH)... (L5). 


0 


. = 
= 
5 I 


u = 
E u r, 
Br} =: 


vs 
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Da sich M(0, 7) und M(l,)) durch die Ableitungen der Verschiebung u(s, 7) in den Punkten 
-8s=0 und s = / ausdrücken lassen, stellt (1.5) eine Integrodifferentialgleichung für das Gleich- 
_ gewicht des Kreisbogens dar. Der Sinn dieser Gleichung ist klar. Das zweite Glied auf der rechten 
Seite, das sich aus (1.3) ergibt, hängt von der Krümmung der nicht-deformierten Bogenachse ab, 
das dritte und vierte rühren von dem Unterschied in den Randbedingungen her.Wie leicht zu 5. 
sehen ist, verschwinden diese letzten Glieder, wenn die Befestigungsbedingungen an den Enden 
des Bogens und denen des Balkens die gleichen sind. Die Gleichung (1.5) geht dann in eine 
FREDHoLnMsche Integralgleichung 2. Art mit regulärem Kern über. | 


Um Gleichung (1.5) in eine FReDHoLMsche Integralgleichung 2. Art überzuführen, kann man 
auch benutzen, daß sich den Verschiebungen y(s, &), die im Punkte D(£) eine der Wirkung der 
Kraft P entsprechende Singularität haben, willkürliche Verschiebungen y,(s, &) hinzufügen lassen, 
die in s stetig sind und in s stetige Ableitungen bis zur vierten Ordnung einschließlich haben. EEE 
Diese zusätzlichen Verschiebungen entsprechen der Wirkung gewisser zusätzlicher Belastungen TR 
an den Enden des Bogens und einer zusätzlichen, über den Bogen kontinuierlich verteilten Be- F 


ar 


a = 


lastung. Die zusätzlichen Hilfsverschiebungen mögen folgenden Bedingungen genügen: 
y0, &) + yı(0, 8) = 0; ud =0...... (1.6a), 


er 20.) 


(0,9) _ 0 , oyll, £) 
E os BR 


ud, 
os i OR 


nn a 


Das heißt: die neuen Hilfsverschiebungen und die zugehörigen Drehwinkel sollen an den Bogen- 
- enden verschwinden. 
; 


Um diese Bedingungen zu erfüllen, brauchen wir nur zu setzen: 
u else tal ma)arals 2 


Damit wird zunächst (1.6a) genügt. Aus (1.6b) und der Gleichung (1.2) finden wir dann a, und a, 
und danach y,(s, &): 


| N 1— s) (l— 
In  U6+9-2@+sAl.... de. 


- Die zugehörige über den Bogen verteilte radiale Hilfsbelastung q,(s, &) wird gemäß (1.3) folgende 
Form haben: 


(58) = a Ss Rn ar (1.9). 
Offenbar ist 
1 S. 


Diese Unsymmetrie der Funktion q,(s, &) spiegelt Eigenschaften der Biegemomente am Bogen 
mit fest eingespannten Enden wider: diese Biegemomente streben für & > 0 gegen Null, aber nicht 
für s— 0, weil in den Endpunkten A und B des Bogens Reaktionsmomente wirken (Bild 1). 
Für die gesamte radiale Hilfsverschiebung v(s, £) eines Punktes des Bogens gilt: 


STE Sa) en ae rn (L10% 
Die gesamte zusätzliche Hilfsbelastung ist 
Ps Eee Se) ee N REE RR TRIER TRIER (1.11). 


“Wenden wir wie oben den Reziprozitätssatz an, so gelangen wir zu der FREDHoLMschen Integral- 
gleichung 2. Art mit unsymmetrischem Kern 


u(E&,n) = v(n, ©) -/ AORNLEINIG ei es (r.127 


isenschaften der Funktion v(£, n) und des Kernes K(s, 8) läßt sich unmittelbar entnehmen, 
n Funktion u(£&, n), die der erte (1.12) genügt, auch die Randbedingungen, d.h. die 
Befestigungsbedingungen für die Bogenenden, erfüllt. Die Lösung der Gleichung bietet ge 
Schwierigkeiten und hat keine prinzipielle Bedeutung. Lediglich auf eine Besonderheit der sic 
aus dem Reziprozitätssatz ergebenden Integralgleichungen wollen wir unser Augenmerk lenken. 


Wie schon bemerkt, sollen die zusätzlichen Hilfsverschiebungen y,(S, 5) in s stetig sein und 
in s stetige Ableitungen bis zur vierten Ordnung einschließlich haben. Im übrigen aber können 


sie willkürlich sein. Wählen wir für y,(s, &) ein Polynom von s von mindestens siebentem Grade, 
11* 
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so können wir stets folgende acht Bedingungen erfüllen: 


v0,)=vl,)=0; a Ds ; 
v0, & _ rl) _ 9. v(0, &) _ Pol) _ 9 
re Te Os de R 


Danach führt der Reziprozitätssatz zu einer Gleichung vom Typ (1.12), die für alle Befestigungs- 
arten der Bogenenden dieselbe Gestalt hat und, weil ihr alle den verschiedenen Randbedingungen 
entsprechenden Biegeformen des Bogens genügen müssen, keine eindeutige Lösung hat: ‚Offenbar 
müssen hier noch die Bedingungen des dritten FrepHoLmschen Satzes herangezogen werden, die 
sich auf den Fall der nicht-eindeutigen Lösbarkeit der inhomogenen Integralgleichung beziehen. 
Bei der Umwandlung der sich aus dem Reziprozitätssatz ergebenden Integrodifferentialgleichun- 
gen in Integralgleichungen muß daher geprüft werden, ob die gefundene Lösung den Rand- 
bedingungen der Aufgabe genügt, und notfalls die allgemeine Lösung nach der FRepHoLMschen 


Theorie aufgebaut werden. 


$ 2. Integrodifferentialgleichungen für das Gleichgewicht dünner elastischer Schalen 


Das in $1 behandelte Beispiel spiegelt im wesentlichen schon die Art der von uns vorzu- 
schlagenden Methode zur Aufstellung von Integrodifferential- und Integralgleichungen für das 
Gleichgewicht dünner elastischer Schalen wider. Der Achse des Bogens entspricht in der Theorie 
der Schalen die Schalenmittelfläche, der Achse des Balkens als „Bild“ der Schalenmittelfläche 
die Mittelebene einer Platte. Als Koordinaten x!, x? auf der Schalenmittelfläche geben wir iso- 
therme Koordinaten vor [3]. In diesem Fall bestimmt sich das Linienelement auf der Schalenmittel- 
fläche aus der Gleichung 


ds? = Fat, 2) (dd + (dd) . . 2.222222. (2.1). 


F(«!, x?) ist dabei eine skalare Funktion der Punkte M(x!, x?) der Schalenmittelfläche. Gleich- 
zeitig sollen x!, x? die cartesischen Koordinaten für die Punkte der Plattenmittelebene sein. Auf 
diese Weise besteht eine umkehrbar eindeutige Zuordnung zwischen den Punkten der Schalen- 
mittelfläche und denen der Plattenmittelebene, und (2.1) stellt den Zusammenhang zwischen 
den zugehörigen Linienelementen ds (Schale) und ds, (Platte) her. In einem willkürlichen ortho- 
gonalen Koordinatensystem auf der Plattenmittelebene erhalten wir !) 


ds® = F*(M) ds? = F{M) GM) Ö® . ........ (2.22). 


Dabei sind die G;; die Komponenten des Maßtensors auf der Plattenmittelebene. Auf der Schalen- 
mittelfläche lassen sich die entsprechenden Größen g9;; (i, k = 1, 2) durch 


Gi = F?(M) G,,(M), Iik = 0 für === k a (2.2b) 
ausdrücken. 


Es mögen nun die Koordinatenlinien x? = z sowohl bei der Platte als auch bei der Schale 
mit den Normalen der jeweiligen Mittelfläche zusammenfallen. Bilden die Einheitsvektoren ?; 
die Basis des Koordinatensystems, so sind die Komponenten des Maßtensors 9; = O(i + 3), 
933, = 1. Die Koordinaten x’ arithmetisieren die Räume innerhalb der Schale und der Platte. 


Die Dicke 2 h der Schale sei gleich der Dicke der Platte Der Einfachheit halber setzen wir A 
als konstant voraus, obwohl die aufzuzeigende Methode ohne Schwierigkeiten auf Schalen ver- 
änderlicher Dicke erweitert werden kann. Ebenso setzen wir voraus, daß für Schale und Platte 
die elastischen Materialkonstanten übereinstimmen. Weiterhin benutzen wir die KIRCHHOFF- 
Lovzsche Hypothese von den ‚unveränderlichen geraden Normalen‘ zur Mittelfläche. 


Wir kommen nun zur Aufstellung der Integrodifferentialgleichungen. Wir bezeichnen mit 
VailN, M) diejenige Verschiebung des Punktes N der Plattenmittelebene in Richtung der 
Koordinatenlinie i, die von einer im Punkte M angreifenden, parallel zur Koordinatenlinie & 
gerichteten konzentrierten Einheitskraft hervorgerufen wird. Bekanntlich ist 


UN, M) = gelM, N) . 2 2 


Die Funktionen v«y (N, M) betrachten wir als kovariante Kom / i 

ai N, A tre r al: e ponenten des Verschiebungs- 
vektors auf der Schalenmittelfläche. Das Kräftesystem, das die Verschiebungen Vo (N M) in 
der Schale hervorruft, besteht aus einer gewissen eingeprägten, im Punkte M angreifenden Einzel- 


1) Hier und auch weiterhin benutz ir di I e b 
) 3 ‚zen wir die von der Tensorr chnun er be i i 
EL HA rohe g h kannte Bezeichnungsweise 
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kraft, einer eingeprägten, über die Schalenmittelfläche verteilten Belastung und aus den Reak- 


tionen. Die Einzelkraft, die den Verschiebungen v,,; in der Schale entspricht, bestimmen wir 
für «= 1,2 mittels des Hookzschen Gesetzes, für x = 3 mittels der in den Verschiebungen aus- 
gedrückten Gleichgewichtsbedingungen. Auf diese Weise gelangen wir zu folgendem Schluß: 
eine an der Plattenmittelebene angreifende Einzelkraft vom Betrage 1 wird in 
eine auf die Schalenmittelfläche wirkende Einzelkraft abgebildet. Der Vektor der 
abgebildeten Kraft wird näherungsweise durch die kontravarianten Komponenten 


EEG) Reh) 


PN NER EEE TER ER 
i | 0 (em. 
bestimmt. Dabei ist 93; = Gy = 12). 


Die auf die Schalenoberfläche wirkende, den Verschiebungen v(,; entsprechende Be- 
lastung kann auf Grund der Hypothese der sich nicht deformierenden Normalen durch ein auf die 
Schalenmittelfläche wirkendes, statisch äquivalentes System von Kräften und Momenten ersetzt 
werden. 


Mit K&, und H£, seien nun die kontravarianten Komponenten der an der Schalenmittel- 
fläche angebrachten, den Verschiebungen v.,; entsprechenden Kräfte- und Momentenbelastung 
bezeichnet; weiterhin mit Sc‘, und L&yz die kontravarianten Komponenten der denselben Ver- 
schiebungen entsprechenden, aber auf den Rand eines Elementes der Schalenmittelfläche wirken- 
den Kräfte und Momente (s. Bild 3 und 4). Die Richtung der Kräfte und Momente und die Be- 
deutung der verwendeten Bezeichnungen sind aus Bild3 und 4 zu ersehen. Der in Klammern 
gesetzte Index (in Bild 3, 4 fortgelassen) kennzeichnet auch hier die Richtung der eingeprägten 
Einzelkraft, die die Hilfsverschiebung hervorruft. Wird der Rand der Mittelfläche nicht von Ko- 
ordinatenlinien begrenzt, so bezeichnen wir die an ihm angreifenden Kräfte und Momente mit S&y 
und L{,. Ihrer physikalischen Bedeutung nach sind diese Kräfte und Momente die den einge- 
prägten Kräften das Gleichgewicht haltenden Reaktionen der Schalenbefestigung. 


Bild 3 


Die den Verschiebungen v.y; entsprechende Belastung wird auf zweierlei Art bestimmt: 


1. Ist = 1,2, so erlaubt es das Hookzsche Gesetz, die Spannungen auf den Grenzflächen 
der Schale zu ermitteln und sie danach auf die Schalenmittelfläche zu übertragen. Zu bedeutenden 
Vereinfachungen führt dabei die Hypothese von KIRCHHOFF-Love& [8]. 

9. Bei x = 3 verwendet man ein elastostatisches Gleichungssystem aus der Schalentheorie 
[4a—c] als Ersatz für die Gleichung (1.3). Das System der Verschiebungen v«y;(M, N) und die 
zugehörige Belastung der Schalenmittelfläche betrachten wir bei der Anwendung des Reziprozi- 
tätssatzes als System von Hilfsverschiebungen und -kräften. 

Nunmehr sei ug. (M, N) die kovariante Komponente des Vektors derjenigen Verschiebung 
des Punktes M der Schalenmittelfläche, die sich unter der Wirkung einer im Punkt N angreifenden 
und in Richtung der Koordinatenlinie i wirkenden eingeprägten Einheitskraft einstellt. Die 
kontravarianten Komponenten dieser Einheitskraft werden durch die Gleichungen 


20 0 (k#i) O3) 
geliefert. 


2) Der genaue Ausdruck enthält außerdem noch Glieder der Ordnung h?k}, wobei k; die Hauptkrümmung 
der Schale ist. Diese Glieder sind bekanntlich gegenüber 1 vernachlässigbar. 


= u. i 
u er a el u g- a 
Die auf den Rand eines Ele 


Lage und die Bedeutung der 


] Momente 
Bezeichnung: 
gelassene Index in Rn. a z 
bungen hervorruft. ich wie oben werden wir die Kräfte und N 
der Basic der Schalenmittelfläche wirken, die nicht mit ur onlir 
mit TE und M& bezeichnen. Die am Rande der Schalenmittelfläche angreifen. 
Momente sind Reaktionen der Schalenbefestigung. Das System der Verschieb R 
der zugehörigen Kräfte betrachten wir bei Anwendung des ei 
Berücksichtigen wir alles Gesagte und wenden den Reziprozitätssatz auf das Haupt- 
Hilfssystem an, so erhalten wir: > 


gs entsprechenden Kräfte und N 


Ugs(M, N) = FM) V&x.(M) KüN) vu, M) a 
— [f FM) V6,,(M) [K%Q, M) uo;@, N) + H3,Q; M) 00;Q Nds 
(8) ! 


+ $ FM) VG,,.(M) [T5, N) vw(Q M) + M6@; N) yQM) . 

— S4,Q, M) un, N) — L4(Q M) oa, Nds nern BAR 
Das Integral ff... dS, wird über die Schalenmittelfläche, das Integral $...ds, über ihre 
3 “ 


(5) 
geschlossene Randlinie erstreckt. 


j Mit M{, und L{, sind die Haupt- bzw. Hilfsbiegemomente bezeichnet worden, die auf den 
Rand der Schalenmittelfläche wirken. Die von den Drehmomenten herrührenden Glieder sind in 
den Hauptkräften 7, und den Hilfskräften S/, enthalten [5]. ®4; und Yo; Sind die zum 
Haupt- bzw. Hilfssystem von Verschiebungen gehörenden Drehwinkel der Hauptnormalen der 
Schalenmittelfläche. 


Nehmen wir an, daß auf der Schalenmittelfläche die Koordinatenlinien x mit den Krüm- 
mungslinien zusammenfallen, und benutzen wir die Hypothese von den nichtdeformierbaren 
Normalen, so finden wir: 


w 
u 


Bm. 


Be " [N 


j 
OK) = D (Ozugyg —. Ozliiya) = Ogligyg u k, Una” a rat er (2.78), 
1 
Ocuyag = 6 (Ozugyı 7 OU) = — O,Ugyg — k, Urea % = Er (2.7b). 
Dabei haben wir die auch weiterhin benutzte Bezeichnung 2; für = und X; für die Hauptkrüm- 
i 


mung der Schalenmittelfläche eingeführt. 

Da H%) #0 ist, geht &g); nicht in die Gleichung (2.6) ein. Analog zu (2.7a—b) werden die 
Drehwinkel y,,; bestimmt. 

Die Gleichungen (2.6) bilden ein System von Integrodifferentialgleichungen für das Gleich- 
gewicht auf Grund der linearen Schalentheorie und erlauben die Bestimmung der Komponenten 
Ugya des Greenschen Tensors (der Einflußfunktion). Diese Gleichungen sind das Analogon zur 
Gleichung (1.5). Durch partielle Integration gelingt es, aus dem Integranden des Doppelintegrals 
die Ableitungen der Verschiebungen u,,; zu entfernen. Machen wir dann noch von der Willkür 
bei der Wahl des regulären Teils der Hilfsverschiebungen v(,,; Gebrauch, so können wir auch das 
Linienintegral $ .... ds, vereinfachen und in einigen Fällen sogar zum Verschwinden bringen 
[4a, c]. Eine derartige Umformung der Gleichung (2.6) entspricht der Überführung der Gleichung 
(1.5) in (1.12). Die durch (2.6) definierten Einflußfunktionen ermöglichen es, die Verschiebungen 
zu ermitteln, die von einer beliebigen Kräftebelastung der Schalenmittelfläche verursacht werden. 


Die weiteren Ausführungen über die aufgezeigte Methode machen wir an Hand eines kon- 
kreten Beispiels. 


$ 3. Beispiel. Sphärische Kuppel 
Zu (2.6) analoge Gleichungen sind bereits auf Fälle angewendet worden, in denen 


F(M)=1. 12 2.000 eo 02 


ist. Es wurden die Gleichungen für das Gleichgewicht zylindrischer Schalen [4e—8] und schwach 
geneigter Schalen [9] behandelt. Für schwach geneigte Schalen ist die Bedingung (3.1) näherungs- 
weise erfüllt [1]. Als einfachstes Beispiel eines Falles mit F(M) = 1 wollen wir die Gleichgewichts- 
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ial- und Integralgleichungen für da 


I Dr F > + hi ‘ 
hungen für eine sphärische Kuppel vom Radius R betrachten. An Hand dieses Beispiels 
lassen sich gleich gewisse allgemeine Besonderheiten der Methode erläuten. 
0 Wirlegen den Ursprung O eines cartesischen Koordinatensystems in einen der beiden Pole pi 
der Kugel und lassen die Achse Oz mit dem zugehörigen Durchmesser der Kugel zusammen- “ 
fallen. Es seiQ die Tangentialebene der Kugel im Punkte N(0, 0,2 R). Dann läßt sich die stereo- ; 
graphische Projektion der Kugeloberfläche auf diese Ebene nach Einführung der sphärischen 
Koordinaten (9, 9) durch folgende Parameterdarstellung beschreiben: u ee 


0 z=Rsindoosg, y=Rsindsing, z= R+ReosB=2 Ron S (3.2), Bi 


rn: = IR1gS 009; Da) RigS sinQ, HUD IE EA ee #4 (3,3). 
Hierbei ist M(x, y, z) ein Punkt der Kugeloberfläche, M *(c*, y*, z*) seine stereographische Pro- 
jektion auf die Ebene Q. Aus (3.2) und (3.3) ist zu erkennen, daß das Teilgebiet 0 < 9 < 9, der 
Kugelfläche auf einen Kreis vom Radius 2 Rtg 2 abgebildet wird. Das System der Hilfsver- 


schiebungen werden die Verschiebungen der Kreisplatte bilden. 


BED 


Wir betrachten die Quadrate der Linienelemente ds und ds, auf der Kugeloberfläche und 
auf ihrer sterographischen Projektion. Es ist 


Ne 
ds? = Re (de + sin ddp); de — ai 2... Bi. 


cost — 
2 


Daraus folgt 


Danach erweisen sich x* und y* als isotherme Koordinaten. Wir werden aber nicht diese Ko- 
ordinaten x* und y* benutzen, sondern bei 2!= © und &?=9 bleiben. (Aus 2.2a) und (3.5) 
erhalten wir 


-4 j 9\-% 
F(M) = cos S z aa (eo 5 ; Ga =G,, =. R?sin?9 (eo 2 (3.6). 


Um nun das System der Hilfsverschiebungen v,,,,; zu Konstruieren, benutzen wir die be- 
kannten Lösungen der Aufgabe, das Gleichgewicht einer unter dem Einfluß einer Einzellast 
stehenden Kreisplatte zu bestimmen [5], [6]. Diese Lösungen genügen den Bedingungen 


Vys(P; M) = Vgya(P; M) a B (& = ale 2) A in (3.7a). 


Dabei ist P ein beliebiger Punkt der Mittelebene der Platte. Auf dem kreisförmigen Rand dieser 
Mittelfläche gelten die Bedingungen 


2030, M) ; 
MT nennen 7b) 


wenn Q ein Randpunkt ist und n die zum Rand normale Richtung angibt. 


’ Wir bedienen uns jetzt der Hinweise, die in $2 zur Ermittlung der den Verschiebungen v,; 

“entsprechenden Belastung der Schalenmittelfläche gegeben worden sind, und benutzen, um die 
Gleichungen (2.6) für «= 1, 2zu vereinfachen, eine Variante der Kırcunorr-Lovgschen Hypothese, 
wonach sich ein prismatisches Element der Schale mit zur Schalenfläche normal gerichteten 
Erzeugenden im ebenen Spannungszustand befindet.?). Dann ist 


Ku=0; Hi,=0 BEL an But (3.8). 


Die Kerne K{, finden wir — ähnlich wie den Kern der Gleichung (1.5) — aus dem elastostatischen 
Gleichungssystem der Schalentheorie. Die Kerne H/, setzen wir auf Grund der KIRCHHOFF-LoOVE- 


3 ü de eine andere Variante der KIRCHHoOFF-Loveschen Hypothese verwendet [8], [9]. Bei der 
ee lettiypothere machen wir davon Gebrauch, daß die Verteilung der Hilfsverschiebungen in 
Dickenrichtung willkürlich ist [4a], und wählen das System der Verschiebungen V(x)j 50, daß es den Forderun- 
gen der Hypothese exakt genügt. Dabei treten im Ausdruck für die Komponenten der abgebildeten Einzellast 


Glieder der Ordnung h?k? auf; diese vernachlässigen wir jedoch. 
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e- a nd längs 
«(Q N) =0; 
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uga(M, N) = R sin®-1 dcs X | r er 


x (ran (N,M) + Rsin ®, fi T4,(Q, N) DO: M) + M4@, N) Yl0, | 


1 ee 
a} (3.10b). Br | 
(i =1,2,3; $=l, 2) s , ö | ud u | 
Hierbei ist Yyr die Projektion der Verschiebung d,., auf die Richtung der Koordinatenlinie i. 


. . . . . :,.1 . für - 
Die Gleichungen (3.10a) und (3.10b) bilden ein System von Integrodifferentialgleichungen für © 
das N der sphärischen Kuppel. Dieses System entspricht der Gleichung (1.5). Seine 


Lösung kann auf die Lösung einer FrepHornschen Integralgleichung 2. Art zurückgeführt werden. 


02(04,) = cos [od N [F KbP, 2) uni Nas) 


nn 


. Beraten SU machen. Bemerken wir noch, daß die Funktionen ay,(0,24) Zr ale Werte 9 
| regulär sind, so können wir an Stelle der Verschiebungen v\,,; die Verschiebungen 
E Veayy(P; M) = umP,M)—uaQM) » 0: 6.11) 
3 in die Gleichungen (3.10a) einführen und finden 
3 ua(M, N) = Rsin*-!9 cos 2 Vorn{N,M) zu ae Zee 
z 


Es ist unmittelbar zu sehen, daß die Funktionen (3.12) die Randbedingungen (3.9) befriedigen. 
Die Gleichung (3.10b) nimmt auf Grund der Gleichungen (3.12) die Gestalt 


Ugy3(M, N) = Duüs(M, N) — JH Ka(P: M) Uns(P, M) dS> ur AIR (3.13) 


an. Damit ist die Aufgabe für die sphärische Kuppel auf die Lösung der Integralgleichung (3.12) 
zurückgeführt. Diese kann mit bekannten Methoden gefunden werden und wird hier nicht be- 
trachtet. 


$ 4. Fragen der Äquivalenz und der Eindeutigkeit 


Die Anwendbarkeit des Reziprozitätssatzes in der Schalentheorie und seine Verträglichkeit 
mit der Hypothese von Kırcauorr-LovE sind streng bewiesen worden [2]. Daher liegt die 
Genauigkeit der oben erhaltenen Resultate in den allgemeinen Gültigkeitsgrenzen der Gleich- 
gewichtstheorie dünner Schalen. In diesem Sinne sind die Gleichungen (2.6) und die aus ihnen 
gezogenen Folgerungen mit den auf der KırcHHorr-Lovgschen Hypothese beruhenden Differen- 
tialgleichungen für das Gleichgewicht der Schalen äquivalent. 


Die Integrodifferentialgleichungen (2.6) haben keine eindeutige Lösung. Die aus ihnen 
folgenden Integralgleichungen können eine eindeutige Lösung haben, und diese wird dann auch 
die gesuchte sein. Jedoch brauchen, wie schon in $1 bemerkt wurde, die aus Integrodifferential- 
gleichungen hergeleiteten Integralgleichungen für das Gleichgewicht gemäß dem dritten FREDHOLN- 
schen Satz keine eindeutige Lösung zu haben. Daher muß man die nach der aufgezeigten Methode 
gewonnenen Integralgleichungen einer zusätzlichen Untersuchung unterwerfen oder eine Lösungs- 
methode anwenden, welche die Erfüllung der Randbedingungen von vornherein gewährleistet. 


Wie schließlich noch bemerkt werden kann, erlaubt es die angegebene Methode, Randwert- 
aufgaben der Schalentheorie auch numerisch zu lösen. 
{ 


Literatur 
[1] Biracog, B. 3., O6ması Teopnsı 060J104eR, T'octexusnar, 1948 r. 


[2] Tonpuenseüsep, A. JI., Teopua yıpyrux TOHkux o6010YerR, Toctexuayart, 1953 r. 

[3] Raran, B. ®., OCHOBbI TeOoPHH UOBePXHOcTeä, y. I, Tocrexusnart, 1947 Tr. 

[4] Kunpuescknü, H.A. (Kinpyveschknü, M. O.) a) OcHoBHi piBHAHHA Teopii 00010H0R i 
Nenkl MeTonm ix IHTeTPyBauHsn, BÖHPHUK Ipanb iHCTUTyTy MaTtemaTuku AH YPCP, M4,5,6, 


5 be 


L - s r } - MY Bi 1 . 
RR N. 5 De ee a u 
ya av V y Temper: Infinite Multilayer Pate — 


— —— e & 


Ei = Fra) L E i = F a 


‚lu 1940 r. 8) HekoTopsie MeToNbI HHTeTPHPOBAHHH ypaBHeHnü PaBHoBecHn yIPyTuUx O060104eR, 
IIlpnksanHaa marem. u mex.,T. IY, Bsim. 2, 1940 r. c) Hadım;keni meronu BusHayenHn IepeMi- 
IIeHb B IMSIHNPMHYHUX 000JI0HKaX, BÖ6IPHHR IPAuB iH-Ty marem. AH YPCP, % 8,1946 r. 

[5] JIag, A., Maremarnyecran Teopua yupyrocru, OHTNW, 1935 r. \ 

[6] Myexenumsnau, H. M., HekxoTopsle OCHOBHbIe 3anayı MaTeMarnyeckoü Teopun ynpyroctu, 
AH CCCP, 1949 r. f 
[7] IIpusanos, M. M., UnterpasepHtie ypasnenus, OHTNU, 1937 r. 


[8] Pemusosa, H. M., Busnayennn mepemimeHb B UHSIHNPMYHUX 00010HKAaX MEeTOAOM iHTeTpanb- 


HHX PiBHAHB, JAH YPCP, N 3, 1958 p. 
[9] ®panaun, B. H., Haxnosernü, C. M., IIpo crırananasa interponuihepenniabHux PiBHAHL 
PiBHOBATH IOJI0THX 00010H0RK, JIAH YPCP, N: 4, 1958 p. k 
Manuskripteingang: 24. 11. 1958 | 


Anschrift: Prof. N. A. KıLTScHEwsK1J, Ukrainische Akademie der Wissenschaften, Kiew, UdSSR 


- ZAMM 40 (1960) Heft 4, Seite 161—165 
Steady Temperature in an Infinite Multilayer Plate 
By VAcLav VoDIORA 


Die Arbeit gibt eine wesentliche Verallgemeinerung früherer Ergebnisse des Autors [1]. Abgesehen von ihrem 
unmittelbaren technischen Interesse zeigt die Arbeit einen schönen Weg zur Bestimmung stationärer Tempera- 
turfelder in langen mehrschichtigen Wärmeleitern. Diese Fragen werden in einem weiteren Aufsatz behandelt 
werden. Die gegenwärtige Arbeit besteht aus drei Hauptteilen: der vollständigen Lösung der beiden Funda- 
mentalfälle sowie der Angabe einer geschlossenen Lösungsformel für den allgemeinsten Fall. 


This paper gives a substantialgeneralisation ofthe author’s preceding results [1]. Inaddition to itsimmediate 
technical interest the following article presents a fine approach to the determination of steady temperature 
fields in long multilayer conduits. These questions will be treated in a separate paper. The article contains 
three main sections: the complete solution of two basic cases and the end closed formulae for the most general 
F problem. 


B paöoTe naeTcA 3HAYHTEJIBHOE 0000IMeHHE Pe3yJIbTaToB [1], paHee MOJIyYeHHbIX ABTOPOM. 

IIomuMo HellocpelCcTBeHHOTO TEXHHYECKOTO HHTEPeca pa60Ta MaeT CTPoÜHbBIaH MeTon MIA 

| oNpeeJIeHuA YCTAHOBUBIUNHXCH NOJIeH TeMNEPATyPbI B CIIONCTEIX TEIMIOIPOBONHNKAX. ONHARO 

| 3TH BOIPOCHI ÖyAyT PaccMaTpHBaTbca B NanbHeümeii pa6ote. Hacrorman pa60Ta Comep;kKuT 

cAIeNylIMMe TPM TIIABHBIe YaCTH: TOJIHOe PelleHme O00HX OCHOBHEIX CIIyyaeB WM BbIBON 
hopMyJIbI B KOHEYHOM BuNeE IA HAN6OJTee OÖIMero cIyyan. 


$1. The first basie problem and its solution 
1) Mathematical statement. Let 
Denen ehren. > 0 0 ee 


be given constants. Let the function f(x) with the period 2% be defined for all real values of x 
and developable into a Fourier series in the interval O<x<2r. Then we seek the solution 
vs =u(&y); i=1,2,...n of the following problem: 


au au — oo <t<+m, Y<y<yiarı; a RER A Et > A N 
Gr. Soy® 
us) =w(c +2 Yy, -e<ı<to,y<y<yys; i=12..un. (2) 
ee oe ee ch ee a ale, 
GE + haea (m u) = + Al uU) = 0, << Hm, 

ydYaı; a 
ee ee ee en 
Physical sense of these relations is obvious from the elementary theory of heat conduction 

in solids. 


2) Particular solutions. Writing the condition (3.1) in the form 


eu) IE) = 2.(a, 05, Kur. ba.) ee. EN 2 on (&)5 
2n 27 
1 a 1 cos k 
Ay = 5n | Ko) do, | Ko) un. dor Bean u 2 lat) 
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states the grounds of assuming the expressions | 
Dr = Xu) Ya) eh na a a US, 
X=l Yo=Cuoy+Do» I=1h2%. „nl; Ag: a en (5.1), 
X; = Ay,coskx + B,sinkz, Yin = Guck km) + Dashk— — 1) 
i=1,23,..,n—l,k=1,2,. 2725. 
Xyr = Anrcoskz + Bu,sinkrx, Ynr=Sshk(yn4ı —Y), a er 


for particular integrals of the set (1); A;. Bj» Cu» Di, are constants of integration. 
The expressions (5) satisfy the conditions (1), (2) and (3.3); they will Br ala 2) if 
and only if 


hirı Kira) Yirlyirı) = hizı Xır1,.(@) Yir,.(Yirı) 
7 Yan) _ Fk Yin), i=1,23,.%,n—- 1, k20 12 Ze 
Yıryırı) Yirı,.(Yi+1) 
where primes on Y mean ordinary y derivatives. 
From this we obtain, using also (5.1)—(5.3), in the first place 


en (ven ” 2) =—l+higı = + A - 


0 
ber) 2 ae (7.1), 
, De ve 
hr On ze In nos hn (v. + a) =—1-— hy (Yazı ie: Yn) a ER a (7.2), 
and secondly 
h; CG;shk(yyı —y) + Di, chk yisı — Y; 
EN je ‚ Cir sh k (Yirı De ‚ch k Yyirı — Y;) Xu@) 
i+1 i+l,k 

i=1,2,. ale ad re 

ln Caı1,25hk (a — Yn- see — Un 

Ray  En=ı, x Sh k (Yn—Yn-ı) + Da-ı,z Ch K Yn— uni) x Ken 


hr ch k (Yn+ı a Y„) 
kK=1,2. 0 RE 


hrı, CauchkyHn—Yy)+DasbkYyaı—yY) _ hir Ci+1,x 


k Casbkywm—y)+Dachkywa—y) x k Dass Pi 
al, A ek 2 TER 
hy no, »chk (Yn—Yn—ı) + D, 1.:5h k(yn— Yn-i) | h,. sh rang —Yn) 
k KR k sh k ( m) + DT k chk(y,— 2 a 2 ; rn k chk(yn+1ı — Yn) ; 


KT, 2 
Now, (7.1) and (7.2) give N% 


Co=H#H: GC: an 


” ka, wc he Se A 
Pe , i > ae he » n—1]), M= (8.1) 
on the one hand and 
Do = Gr Co = Ho Cho: = 1,2 ee 


w=g| 2 2* ‚Ale ıQjo + Ha Ynı ), i=1,2,...,n—2; In On + 7 Yarıs 
adE 


i\j=i4 


h 
Qu =— +1 1)y. Tann Are 


on the other. Similarly, we easily get from (7.5) and (7.6) 
Da = 440} i=12,..,2—1,k=12. ya (8 
_ hir ‚th k(yi- +1 —yı) —di+l,k Ih: ++ kthk (Yizı —l, 


Ye 
hir — 9i41,» [K + Auyıth k Gas Hy] 
i=1,2,...,n—2,k=1,2. 
ee ur kuthgn— Yn)]th.k u rd 


k ar An th k 7 Zi Yn--ı) + h}, th k (Ynzı — Be k= 1, 2, ... (8.3.1) 
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and, after that, the equations (7.3), (7.4) leäd to 
hirı Cr 


Kr. @) = ae [sh k Yi4ı —Y) + Gr ch k(yi4ı — Y)] X) 
Il, 2,......0n—2, k=12.., 
er Dr Cn-1,% - 
Kr) = — h„chk (yn.ı— yn) [sh k (ya — Yn-1) + m-1,2Ch k (yn — Yn—1)] Kn—,2@) ; 
IE RATE 
From this we deduce without difficulty 
| BG 
Re N ne ee 
M:+1,2 Ci+1,% 
H2G 
a 1k = 
nk(X) M,_ x chk (na —yn) Nr Kıı(®) » k=12,..., 
i+1 
Mita I 9; De RD ne Bell) 2, Br 


Nir1,k = an k (+1 En y,) + Apr ch k (y+1 = y,)] ’ 
= 152,5 2, a „ee 2a EB 
Using (5)— (5.3), (8.1)—(8.4) we finally get the following system of particular solutions 
Dig = Virl@, Y): 
Io = H; Co (y 7 io); = 1,2,..., n—1, Uno = Hn Cyo (YUn+ı — Y) ET (9.1), 
Djr = Pix Cır (A1r c0skx + B,,sinka) [(hkYy— y) + gxshky—y)], 
Deren Ken oer, 


OA er Binde kl, 2%: 


H;N; EN 
er adeoaee ne NW era I 
a. at M;: ( ) Eat MM, 5cCh K (yy41 U.) 
k=1,2. 2202 m 


The expressions (9.1), (9.2) satisfy all conditions of our initial problem with the exception 
of (3.1). 
3) Solving of the beginning problem. If we assume, as usual, the wanted solution 
u; = U;(x, y) to be composed from the expressions (9.1), (9.2) in the form 
Ur, et BIT Re RT 1 
k=0 


our last unsatisfied condition (3.1) will be expressed — see also (4) — by the relation 
I (a, cos kx + b, sin ka) = 3 v4. Yı) = Co Yı + Mo) + X Cr (Arcoskx + B,sinka). 
k=0 k=0 k=1 

From this we have 


a 
CH == nd = rd A br Kr 1% 2,80% ee (11) 


Yı + 90 
and the solution of our first fundamental problem is, by (9.1), (9.2) and (10), 


ae) + 3 Pix (9 cos ka + bu sin ka) [chk(y—y) + gu sh ky— yo], 


= 1 
ee ee 
BA) =. I, u +2 Par (a, coskz +b,sinkz)shk (min —Y) -»:.:... (12). 
Yı + Qio Pr 


The constants a,, b, are determined by (4.1), the quantities H, by (8.1), the expressions g;; can 
be found from (8.2), (8.3.1) and the values p;; from (9.2). 
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$ 2. The second fundamental problem and its solution 


We now discuss the problem differring from (1)—@.3) merely in the boundary conditions 
(3.1) and (3.3). In their place we require now that 


u y)=0, -o <2e<tm0: vo 2... 000 (13.1), 
1,6% Y) =Fa), -o <<... ne ui (13.2); 


the function F(x) is submitted to analogous assumptions as formerly f(x). 


There is obviously no substantial difference between this new and the previous-problem 
and we easily can obtain the solution of our present case almost without any calculations directly 
from (12). In fact, to do it we have only to replace former quantities by the new ones according 
to the following readily comprehensible scheme: 


U;; Yis h;, hi, ; . 
a Eh ae ee 14). 
at Un, und — In+2-, Fi) (14) 
In this way we find the following solution of our second problem 
A —y 2 . 
1, y) = — = Fr a 2 (A,coskx + B,sinkz)shk(y—y,), 
A,H; ‚Y+tro 


u, y) = — 


TE + 2 bı(dı coskx + B,sinkx) [ch k (y;}ı — y) 


— Tr, sh k (yiyı —Y)] ; i=23,.:., oe 
Y-+ Tno 


Bee + tr 
n-- no 


+2 (A,coskx + B, sin kzx) [ch k (y„..ı—Yy) — Ta; Sh k (yn+1—Y)] 


(15). 
The symbols and quantities occurring in this result are determined by the following formulae: 


n 


2 an 

A 1 coskw 

A = ER e NER. N 

0 5 | F@) ao, Be | Fo) Sn dv, k=1,2.2u » 5 E57 
0 0 


i ! r 1 n—1 
Ba (2.8) ‚ 1=2,3,...,n—l, Pao= zp (Hn Ro— 3 H, Ro +), 
n »—=2 


1 h; 3 
Ro= + 1) ER N, © (15.2), 
„aut lktbthkg—yltbkg—y) Eee 
k+hthk,—y) + Mk, thk—y,) ’ Er 


ne h,th k(yyı —Yı) + Ta, [ki + kthk (yıyı — y)] 
+ Te [k+ kithk(yıyı — Yı)] er 


3,4 Akira 


N, 
es — a ee BE 
Hn My,» ch k (9, — yı) 
H; NaH-i k 
I TG ET EDS Bus en 
H,My41—,k u. 1, ke 
m n—l 
M, = I, T; ’ j =D . a 
ik ae jk I 2,3,...,n—1l,k=1,2, 


N;r — (— N li [sh k (y;4 ı—Y)—Lyr ch k(yy+1—y,)] > i=2, 3, .s., A, k= 1; 2, en. 
(15.4). 
$ 3. General case 
Our last problem contains both pr i i i 
a ae re cases as the special ones. It differs from (1)—(3.3) 


superposition of the results (12), (15). now by (13.2), and its solution follows evidently by 


# 
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We then have 


ar Y-+ go S : 
U(%, y) = qy eg +2 (a; coskx + b,sinka) [hky— y)+q,shky— yı)] 
Mal en tu (A, coskx + B,sin ka) shk(y— y,) 
Ich, Un-+1 + Tn0 k=1 a 
3 Pi H:; Yy % Gio S . 
u(& y) =@aH; —— +3 Pir (a, coskx + b,sinka) [chk(y— y,) 
Y+ MW a1 


AH; _ Ytro 


+ g,shky— y,)] — H IRRE ER 
N N no 


+2 b;; (A, coskx-+ B, sin kx) [ch k(yyı —y) —r;, sh k (yy41ı —Y)], 
z Den 


Unit, 0) = 0,H, au +2 Pnr (9, coskx + b,sinkx)shk (ynrı —Y) 


YTtQo = 
A F y + Tn0 = A . 
+49 — +2 (A; coskx + B,sin kx) [chk (y„n41 — Y) — 25h k (yn+1— Y)] 
a Yn+ıt no 1 


(16). 


This is the complete solution of our general problem, expressed mathematically by the relations 
(1)—(@.2) and (13.2). 


For the sake of better clearness and applicability we again remind of the signification of 
all symbols and abbreviations occurring in (16): a,, b, are defined by (4.1), A,, B, by (15.1), 
the values q;, can be found from (8.2) and (8,3.1), r;, from (15.2) and (15.3), p;, are determined 
by (9.2) and Z;, by (15.4). The auxiliary quantities H,; are given in (8.1) and the expressions 
M;% N;; in (8.4). And finally, ch, sh, th mean hyperbolic functions. 
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Stationary Temperature Fields in Multilayer Cylindrical Tubes 
By VAcLAv VODIÖRA 


Ziel der Arbeit ist die Verallgemeinerung der in [1] angegebenen Resultate für den Fall eines unendlich 
langen hohlen Zylinders bestehend aus n Schichten. Die Arbeit zerfällt in zwei Teile: zuerst wird der Fall 
des Kreiszylinders und dann der des elliptischen Zylinders betrachtet. 


The object of this paper is to generalise the results given in [1] for the case of hollow infinite cylinders 
composed of n layers. Our considerations are divided into two parts. First comes the case of a eircular and 


next that of an elliptical eylinder. 


3anaya pa6oTbI COCTOHT B OÖ0ÖIMeHHN PesyJIbTaToB, MAHHBIX B pa6orte [1], ma cıyuaä 
ÖecKOHEYHO MNIHHHOTO IOAIOTO IMAMHAPA, COCTOAIMerO U3 n CAOeB. Pa60Ta COCTONUT U3 NBYX 
yacreü: B IIePBofi yacTu paccMaTpuBaertca cıyyali KPYTOBOTO, BO BTOPOH — YIITUITHYCCKOTO 


NHJIUHNPA. 


1. An Infinite Cireular Cylinder 


1.1. Statement of the Problem. Assume that , >>> @+1>0, h,>0 and 
K>0(u=2,3,...,n)are constants and /(d), F (9) functions with the period 2 and develop- 
able in FOURIER series in the interval 0 < 9 < 2. Then we have to find the periodie solution 
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u; = u,(0,9)(i=1,2,...,n) of the following problem: 


271. . : 
en aa er? Oi+1 <o <d5 0zB9 z2n Weßzr dl), 
0 4 


ne, BE), 0a 1.21), 

i d iv v ; 
Sr hir (u, — Ui41) = 2 — kya(y —UH)=0, 

o=an, 059820. (i=12...,n—QL(2Y, 

(1.2.3) 


Un (On-+13 d) = FiD), 0 sas2n Pe a an 


The physical interpretation and technical significance of these equations are evident. 
1.2. Transformation into a Linear Case. The familiar relations — see Ref. [1] 


E=d, n=inD u =. (1.3) 
transform (1.1)—(1.2.3) into the following more simple problem 
u, Au; s 
je + ni =0, oo <i<to, Mm<n<mNa e= 12,2) 2 2 
En) =w(+2n,n), -o<f<too,n<n<pa (i=1,2,..,n) (15), 
WED) je), oc too. (1.6.1), 
ou; r7 o i 7 
en + Hy (u —ui41) = +Hnıw—un)=0, —o<g<+m, 
0. M+ı G=1,2%...,Aa—- We. 020 
ln) =, FREE TR HE A RT (1.6.3), 
where 
7, = In | ’ Hıyı = hiyı0i41> His = hizı Bi+l 27» Pa (1.7). 


1.3. Solution ofthe Problem (1.4)—(1.6.3). There is merely a difference of notation 
between the problem in question and that of Ref. [2]. The correspondence between the two 
notations is as follows: 

% Y, Yı Yu» Bir, Air, | 
5, N» 0, Nis Hrn, H;+ı | 
Thus, the wanted solution u;(&, n,i=12,...,n of (1.4)—(1.6.3) can be immediately 


obtained from the correspondi tions i i i 
ae ponding equations in [2] by simply interchanging the symbols, accord- 


“mw WE ER Eee ae 


u(&,n) = «li + 2) +2 (a cosk&+b,sink&)(chkn + qushkn) 
% A, ri 


H, BE er — Ihe (A, cosk£ + B,sink&)shkn, 


n+M — 
u(&n)= a H, 1 = - +2 Pir (a, cosk& + b, sin k &) [chk(m —n,) 
A H,; Y +rn > 
0” eg + % be (A, cos k& (1.9). 


By innen 
+ B, sin k£) [chk ya —m)— Tr, Shk(mı—n)] (i=2,3,...,n-—1), 


+ gs km —n)] — 


e Nn+1—N = 
u,(E, = rg ee u Seh 
0), “ + 2 Pnx (a,cosk& + b,sink&)shk (Nn+1ı — N) 
N + 170 
le 
Nn+1 + Tno 
— Tnk sh k (Mn 07 n)] 


+ A +23 (A, cosk& + B, sin k&) [ch k (NMn+1— N) 
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en Bi quantities and abbreviations uBeR.DE here are determined as follows — see also 


2n 
r | 
=) Ko) do, nie ): er (Kelle EIN: 


u 


2 27 
=] A 1 k 
N Ende, -1[ ee & 

ar, (o) do, B; r, F(o) sin ku I® (k E32, OR 29,2% 


dio = m, SH il Q,0+Hn +1) (i a n—2), An-1, o—=oro hi — Mr 


ne Hith k miaı — N) — %r1,3 [Air + kthkyrı — m] 
Hi — gie [k + Hizı thk (mai n)] 
Bells ee, nz e2 Ri, 
Hn + [k+ Huth k (+1 — m)] th k (fm — Mn) 


(1.93), 
k+ Hnthk(m—M-ı) + Hy th k (Mu4ı— Nm) ER a 


An—1,k — 
EN h.ht...h, 
ee N a Br nt NER TL IST. 
1 a re a el ae ein 2 


= 


A i h 1 n—1 
Tjo RT. ars, it Le Tun = 5 (Ha Ban — 2,4, Ru); 
 H+(k+B,thkn)thk(g— n,) 


Toy = - = Kr, ; 
7 k+H,thkn + H,thk(m —n,) \ 
Be un Wr: [H: + kthk (na — m] (1.9.4), 
i z LA 
H; + ri, [k+ H;ithk maı — n)] 
(==93,4,. 3 k=1,2, )3 
1 h, 
RBELLEST (H=2,8, ..., Ar), Pnk = HnNnr - | 
2 M;: Mizi:E ch k (Nn+1 — 1) 
(Kerl), 
i+1 a ( e 
Mi = I] 9: een 2.) (1.9.5). 
Nir, = II [sh k (41 —n,) + Gr ch RK 41 — 9,)] 
2 ee BR A ern 
N BEN, 
hr = a]; ‚= _ u Nnrı—in (2, 32.00: k=1,2,2 
I; Ma—1,%k ch kn H NEE ji 
% n—1 
Mu= I,’ TEEN Tl KEN, 2 u \(1.9.6). 
ı=n EV 
Nr = (1)! 112. [sh k (m 41 — 9) —Tju ch k (n+ı — E e- 
= (=2,3,. erlen, 


1.4. Solution of the Problem (1.1)—(1.2.3). This solution is given by the formulae 
(1.9) in which £ is to be replaced by®andnby In (g,/o), according to (1.3). Allthe remaining quan- 
tities and auxiliary expressions are to be caleulated from (1.9.1)— 1.9.6). 


a a nn 5 
u R 


» >% 
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2. An Infinite Elliptieal Cylinder 
9,1. The Statement of the Problem. Let D,(i=1,2,...,n) be the domain between 
the confocal ellipses e,, e;ı which are determined by 


detpp—pg=0, p—g=e, c>0, cconst. WE1,232..90+ 0.00 


Hereat we assume 
0 <Ppı <Pa <''' <Pn <Pati» 
so that we have also 
I<n<pR< <m <Mrı: h 
Let p(x, y) denote a function prescribed on e, and D(z, y) another function prescribed 
on &,.+1. In the system 


z=cceosechn, y=ctsnöshn ..., . zu. ns (2.2) 
of elliptical coordinates £, 7 our ellipses (2.1) are given by 
Er erT, BEN un (=1,2,...,2n+1l) -. .#3) 
Ze De 


We have especially 
p=cchn, g=eshn; z=pcosf, Y=qsmnEe..... (2.4) 
along the curve e, and all subsequent calculations will be performed under the assumption that 
the periodic functions 
Plz, y) = P (pı cosE, sind) = /(&) , Dia, y) = P(Pu+ı 6088, Qu+ı sind) = Fl) . . (25) 


can be expanded into FOURIER series in the intervalO SE <2n. 
fErR,>0,,>0(i=2,3,...,n) are given constants, we are interested in finding the 
solution , = u,(&, y) (Ü=1,2,...,n) ofthe problem (whose practical meaning is quite evident) 


u u, j } 
Fa a ee En 
1.59) = Plan Aalong er , 2 area Hab RE A BE (2.7.1), 


ou; ou;.. r 
Zn + han (y — ww) = än ! + han (y — Ur) = 0 


along ke =, 3%...,,n—) 7.@72, 
u,(%4y) = Pz,y) alone 2 N. „ere (2.7.3). 


The derivative ö/On in (2.7.2) is taken along the normal to e;;,, directed from D; to D;,;.. 


2.2. Transformation into a Linear Case. In order to avoid the well-known mathe- 
matical difficulties which would lead us far into the theory of Mathieu functions, we confine 
ourselves to the types of tubes for which ch? 7, > 1— see also Ref. [1]. Then the Egs. (2.2) 
reduce the above problem (2.6)—(2.7.3) to the following more simple form: 


eu, Au, h 

oe Fe ‚, -o<{i<+to, „<y<na (el; iu cn) (2.8), 
u(&n)=u,(£+27,n), —oo<f<+oo, H<N<NH 1; Den n) (2.9), 
uen)=fd): << N, (2.10.1), 


ou, ou;. f 
F7 + hr PiH (u u) ne + kHhıPprı wu) =0, 
-o<f<+o, yanyı eisen (2.10.2), 
Uy(&, Nn+ı) = F(E), Zei >, <ö<+to a 6. Mu mu. Fr (2.10.3). 


2.3. Solution of the Problem (2.8)—(2.10.3). The only di£ 5 
from the case treated in Ref. [2] is one of Aninhe table a 


yy, bh, ki 
nn, hp, K Pi 
illustrates the correspondence between the two notations. 


oA 


v 
% 
n 
ü 
n 
» 
# 
2 


yr a ee 5 
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 Applying the substitution (2.11) to the resulting formulae of Ref. [2] leads directly to 
the following solution of the present problem: 


EU S 
u&n) = eg rn (a, cosk& + b,sink&) [(hkm—n) + Qrshk Y— nn) 
A, nom 


FE m er Ar soskr } Bein köychkg m), 


u, 7) = aM; " ee + I Prrlay cos k& + by sin k&) [ch ky —n) 
Q10 Bi 


AH, n-+ro 
ıg sh k =) Or Ü 
nn, (2.12). 


+2 br (Ar cosk& + By sin k&) [ch ku —9) — rue sh km —m)] 
(Be 28 en 


u„(& 9) = % Hn .- Sr > Pnx (4; cos k& + b, sin k&) sh k (m +1 — N) u 


+2 (A, cosk& + B, sin k£) [ch k (m+ı — N) — Inz Sh k (m+ı — n)] 


Allthe auxiliary quantities appearing in (2.12) can be calculated as follows — see again 
Reis [2]: 


sinko 


2an 2n 
went . | coskw = 
= uf [w) do, = z, ®)- .. do RE Ne 
0 


et Aral | coskw er 
— nr F(o)do, Br F(o) le do (Kl 2 122 
0 0) 


1 2 . 
dio - | 2 Hm Qi + Hn Ma) (=1,2,...,n—2), 9m-1,0 = = ty mal 
IE 


hirı Pizı th k (Mi — N) — 4ir1,# [hir Piyı + kthk (mr — m)] 


RE hizı Piz — Qi+1,x [K + hizı Pizı th k Mi — m)] 
Vai sin ea 1, 2 
_ Mn Pa + Ik + fu path k u — lt km—m)  K-1,2,..,, 12123), 
Bet k+ An Pn th k (Mm — Mn) + Fa Pn th k (Mn41ı — 9m) 
Bahr, Di er re BEE a 
RR N TUR, (>02, Ba ass 15 = Be EL 
= E ni 1), ER) 
hp; \ı 
ki > nt) 
T;o Ser: DDR, Nı (= 9 
i w>2 
1 n—l ) 
Tnoe = H. (HM, 2 ZH, R,+ m) » 
& = 
h, Pa A [k + h; pgth k (ns — ı)] th k (ns — 93) (k ae. 1: 5 er 3 
KL RPpthkm—n) + Kpethk (m — N) (2.12.4), 
BR Pi thk ner 1) + R-ı,k [h; Pi + kth k (Ni+ı u 
1 hp -+r-ıelk+ hip th k (m+ı —m)] 
insg 
1 (& i 
. IE 11% WE Do. Han 
R;o h; m h; N ( ) 


170 
H,„N 
No) ent 0 ie EB ee 
Dik —= ne (i = 2, 3 en Lin 1) ’ 17 M„_1,chk (Mm +1 — N) ( 
? 
‘+1 
PR i=1,3,..,0n— 23; k=12...); 2.12.5), 
Mi+1,% II%: ( i 
Narr = IT Ih ka — m) + ech k mr —m)] } 
+1,k II ((=1,23,...,n—l;k=1,2,...) 
t = Nat ’ En (i=2,3,...,n—1; k=1,2,...), 
& 0:2 M„—ı,» ch k (N — N) H„Ma+ı-i,k 
Wie Tem keaseeeınıe 
I | (212.6). 
Na = [sh k (n;41 — 9) — Tjr Ch k (m — M)] 
a (= 23. hm kohle 
Ip a 1 
=] oe ie 1 a N ben ) 
“2 "4% 


i ired i laced by x 
lution required is expressed by the formulae (2.12) where & and n should be rep 
El y, ii to 23). Of course, as we have seen, our solution (2.12)—(2.12.7) is merely an 
approximation, the degree of exactness increasing simultaneously with 7,. 


found quite similarly as in Ref. [1]. 


111 v. 


[2] V. 
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gleitende Rad wird an der Stelle der Berührung mit der Schiene abgeschliffen, und es entsteht 
an ihm eine sogenannte Flachstelle. Bei weiterer Rotation verursacht dann ein solches Rad eine 


stavebnich konstrukci, III. &äst, Vybran& state“ (Baudynamik, III. Teil, Ausgewählte Aufsätze), welches im 
Staatsverlag d. techn. Liter, in Prag i. J. 1960 ( 


L. FrysA, Schwingungen des unendlichen, federnd gebetteten Balkens ... 


fd = p(pıecose, sind), F(&) = D (Pu+ı 60SE, Qn+ıSinE) .. . (2.12.7), 
9.4. Solution ofthe Problem (2.6)—(2.7.3). Assuming that n, is sufficiently great, the 


Another approximate solution, which can also be of use for practical purposes, can be 
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Schwingungen des unendlichen, federnd gebetteten Balkens 


unter der Wirkung eines unrunden Rades”) 
Von L. FrYBA 


Der Einfluß einer Unregelmäßigkeit am Rad oder an der Fahrbahn und besonders der Einfluß des Rades 
mit einer Flachstelle auf die Schiene wird untersucht. Das Problem wird als erzwungene Schwingung eines 
unendlichen Balkens auf elastischer Unterlage und eines Systems von 2 Freiheitsgraden dargestellt. Das 
System von Differentialgleichungen wird durch die Integraltransformationsmethode gelöst. Laplace-Carson- 
sche Bilder werden in Volterrasche Integralgleichungen geführt, die numerisch gelöst werden. 

The influence on the rail of irregularities at the wheel or rail, in particular that of a flat spot on a wheel, 
is investigated. T’he problem is presented as that of forced vibrations of an infinite beam on elastie foundation 
and of a system of two degrees of freedom. For solving of the differential equations the method of integral 
transformations is used. T’he Laplace-Carson transforms are transformed into integral equations of the 
Volterra type which are being solved numerically. 

B pa6ore n3yuaerca BIIMSIHHe HEePABHOMEPHOCTH Ko,leca MIU MOCTOBOH MH CIHENHANBHO 
BJINAHNE Koleca C BBIÖONHON Ha PeIpc. Zanaya IpencTaBılneTcH KaK BbIHYIKICHHOE Konebanne 
ÖEecKOHEYHO JUIHHHOH OAaJIKU HA YIPYTOM OCHOBAHUN MH CHCTEMEI € ABYMA CTeIIEHAMH CBOÖOAEI. 
Cuerema nuphepenmmanpHubIx ypaBHenumi pemmaercan mMpM IOMOIIM Merona HHTETPANbHBIX 


npeoöpasogannä. Vlao6parkenns Jlamsaca-Kapcona IIePeBOAATCA B HHTETPAJIBHBIe YPaBHeHngA 
Bonteppa, KOTOPbIE B CBOIO OUePe]B PemmaloTcHh UMCHEHHO, 


1. Einleitung: Geschichte des Problems 
Die unrunden Räder entstehen beim Bremsen der Eisenbahnfahrzeuge. Das festgebremste 


*) Ein Teil aus der Dissertationsarbeit [1], die fast vollständig im Buch von Prof. KotLousek „Dynamika 


tschechisch) erscheinen wird, veröffentlicht werden wird. 


EA Are ni ie WE SE Be Ai be re 


— 
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Zusatzbeanspruchung der Wagen sowie der Schienen. In diesem Aufsatz wird nur der Einfluß 
der unrunden Räder auf die Schienen untersucht. Ihre Zusatzbeanspruchung durch die Flach- 
stelle ist so groß, daß sie Schienenbrüche verursachen kann. Es ist ein Fall bekannt, daß eine 
Flachstelle von 30 cm Länge auf der 4 km langen Strecke 19 Schienenbrüche verursachte (s. [2); 
in der Deutschen Bauzeitung, 1893, S. 155, wird sogar von 81 Brüchen berichtet. — Die Flach- 
stellenhöhe erreicht bis 20 mm, ja sogar 31 mm (s. [2], [7]). Das Vorkommen von Flachstellen ist 
häufig. Rusın führt in [7]an, daß 1,56% von beobachteten Wagenrädern Flachstellen aufwiesen. 
Das bedeutet, daß jeder Zug ein bis zwei Achsen mit Flachstellen besitzt. Aus diesen Tatsachen 
ist ersichtlich, von welch großer praktischer Bedeutung die Lösung dieses Problems ist. 


SALLER [2] (1921) befaßte sich als erster mit diesem Problem. Er setzt die Flachstelle in der 
Form einer Sehne voraus und berechnet die kinetische Energie vor dem Stoß. TIMOSHENKO 
[31 (1926) zieht die Schienenunebenheit in Betracht und vernachlässigt den Einfluß der Schienen- 
masse sowie der abgefederten Masse auf die Schwingung der ungefederten Masse. Die Voraus- 
setzungen im Beitrag von Anrer [4] (1935) und im Aufsatz [5] (1951) sind wesentlich vereinfacht. 
Popp [6] (1952) geht von den gleichen Annahmen aus wie SALLER [2] und berechnet die Biege- 
spannung. Nach Popp wächst die Schienenbeanspruchung mit zunehmender Fahrgeschwindig- 
keit. Dies widerspricht jedoch den Messungen, s. [11]. Rusın [7] (1952) stellt eine empirische 
Formel für die Schienenspannung auf. GoR’Kov [8] (1955) nimmt die Flachstelle nur auf einem 
Rad der Wagenachse an und gibt die durchschnittlichen Stoßkraftwerte an. In einer früheren 
Arbeit des Verfassers [9] (1956) wird nur die erzwungene Schienenschwingung in Betracht ge- 
zogen. — Die experimentelle Erforschung des Einflusses eines unrunden Rades litt zuerst unter 
den unvollkommenen Meßgeräten (die amerikanischen Messungen i. J. 1920, 1937, 1942 und die 
deutschen i. J. 1932 [10]). Umfangreiche Versuche wurden in USA erst i. J. 1947 [11] und in der 
DBR i. J. 1952 [7] durchgeführt. Rusın [12] (1953) verglich dann die Ergebnisse der deutschen 
und amerikanischen Messungen. Auch BerzHoup [13] (1957) befaßt sich mit dieser Frage. 


Aus der Geschichte des Problems folgt, daß sich die experimentelle Erforschung in den 
letzten Jahren dank den modernen Meßmethoden sehr gut entwickelt hat. Es liegt jedoch noch 
keine geeignete Methode für die theoretische Bearbeitung des Problems vor. Die Ansätze in [2] 
bis [9] lösen es nicht vollständig. Die Unmöglichkeit, den zeitlichen Verlauf der Erscheinung 
zu ermitteln (mit Ausnahme [3], [9]), und das Ersetzen der Schiene bei der Idealisierung durch 
ein System von einem Freiheitsgrad (ausgenommen [9]) sind die größten Mängel der bisherigen 
Ansätze. Das Ziel dieser Arbeit besteht also in der theoretischen Ermittlung des Einflusses eines 
unrunden Rades auf die Schiene, 


2. Festlegung des Problems 


Wir legen unseren Untersuchungen folgende Voraussetzungen zugrunde: 

1. Die Schiene wird als unendlicher Balken (von konstantem Querschnitt) auf elastischer 
WINKLesscher Unterlage idealisiert. (Die Brauchbarkeit dieser Theorie wurde in den letzten 
Jahren wieder von BErzHorLD [13] (1957) experimentell bestätigt.) Die BERNOULLI-EULERSche 
Differentialgleichung gilt für die Querschwingung dieses Balkens; dazu nehmen wir Flüssigkeits- 
dämpfung proportional zur Schwingungs- 
geschwindigkeit an. 

2. Der Eisenbahnwagen kann als ein Sy- 
stem mit 2 Freiheitsgraden aufgefaßt werden: 
das Rad mit den übrigen ungefederten Teilen 
(wie Lager u.ä.) bildet den Massenpunkt m, ; der 
Teil des Wagens mit der auf ein Rad entfallen- 
den Belastung ist die abgefederte Masse m,. Die 
Feder zwischen den Massen m, und m, mit der 
Federkonstante C weise wieder Flüssigkeits- 
dämpfung auf. 

3. Der Verfasser zeigte in seiner Arbeit [9], 
daß der Einfluß der Geschwindigkeit c, mit der 
sich der Raddruck auf der Schiene bewegt, für 
die üblichen Eisenbahngeschwindigkeiten ver- 
nachlässigt werden kann. Wir werden also vor- 
aussetzen, daß sich das unrunde Rad mit kon- 
stanter Geschwindigkeit (ohne Reibung) auf der Stelle x = 0 dreht. Wir nehmen dabei an, daß die 
Belastung stets in den Punkt x = 0 übertragen wird. Das ist berechtigt, weil die Flachstellen- 
länge im Verhältnis zur Durchbiegungswellenlänge der Schiene klein ist. 

4. Der Einfluß des unrunden Rades wird so betrachtet, daß sich die Entfernung des Rad- 


schwerpunktes von der Balkenmittellinie um /() ändert. Dies kann auch durch Fahrbahnuneben- 
12% 


Bild 1. Unrundes Rad auf der Schiene 
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. Sg: it 
heiten (z. B. Schienenriffeln) verursacht werden. Die Kurve /(f) wird im Falle des Rades mi 
einer Flachstelle folgenderweise gegeben: | 
ET a et N 202.1 2.08, EEE 
NO für, ES OEH a 


Dabei ist: . 
T die zum Durchfahren der Flachstelle erforderliche Zeitdauer; 
jı() die Flachstellenordinate in der Zeit t (Länge AB in Bild 1), 


“ a=f(T/2) die Flachstellenhöhe, 
l=cT die Flachstellenlänge, werte 
c die konstante Fahrgeschwindigkeit. 


Tafel 1. Flachstellenordinaten 


t(t)/@ | 


re et — Ä 
ELEREIEN Et 
| \ 
| | 0,905 


| 
) | 0,016 | 0,071 | 0,172 | 0,304 | 0,455 | 0,627 | 0,783 


Bild 2. Flachstellenordinaten 


Die Ordinaten /,(£) wurden gemessen. Ihr Mittelwert ist in Tafel 1 und in Bild 2 angeführt. 
Wir drücken die gemessene Funktion f,(£) mit Hilfe einer Fourierreihe als periodische Funktion 
mit der Periode 7 aus: 


oo x 
hd =. —- I! gcsjot= La ll —esjeoh TI a2, 
il je1 


oo 5 1 
wobei; = 1 2, 3 so = 2 =9gellug, =2 4;, weil /(0) = 0; und = T | hd dt, 
j- . . 


v0 


gr 
= — . | fi) cos wtdt die Fourikerkoeffizienten sind. Wir erhalten sie durch numerische 


0 
Integration aus den Werten der Tafel 1: 
a, = 0,4806 a; & = 0,5090 a ; A, = — 0,0202 a ; A, = — 0,0070 a; 


Wenn wir in (2.2) nur das erste Glied berücksichtigen, bekommen wir annähernd 
(a, 8 0,5.a): 


1 = a (1 —e0s 0) » zZ all—cosot) De 


Die Kurve (1 — cos w1)/2 wird zum Vergleich in Bild 2 gezeichnet. Die Abweichungen dieser 
Kurve von der gemessenen f,(t) sind sehr gering (meistens im Rahmen der Meßgenauigkeit). Es 
genügt deshalb /,(f) nach Gl. (2.3) darzustellen. 

Wir lassen die Tatsache, daß sich der Stoß periodisch wiederholt, außer acht. Alle experi- 
mentellen Messungen [1], [7], [10], [11] zeigten nämlich, daß die durch die Flachstelle verur- 
sachte Schienenbiegungsspannung und Durchbiegung weggedämpft ist, wenn die Flachstelle nach 
einer Umdrehung die Schiene wieder stößt. 

5. Die Herrzsche Gleichung für die relative Annäherung a(f) des Radschwerpunktes und der 


RE 3 ee 1/3 
Balkenmittellinie ist «(f) = x’ P?/3(t), wo x’ — 5 2 42 e) rent! der aus den Hertzschen 
DT; 


Formeln [14] für die Berührung zweier unter dem Winkel r/2 gekreuzten Walzen mit den Halb- 


ee EIN NND VUEHE GEN WE, VuE 


FR RE 7: 5 ’ Eu 
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messern T, r, erhaltene Koeffizientist. Die zwei gekreuzten Walzen sind in unserem Falle der Rad- 
reifen und der Schienenkopf (aus gleichem Material, Poıssonsche Konstante v< 1 und Elasti- 
zitätsmodul E). P(t) ist die Berührungskraft. Die Herrzsche nichtlineare Beziehung wendete 
TımosHenko [15] (1913) für seine Lösung des Balkenquerstoßes an. Er erhielt eine nichtlineare 
Integralgleichung, deren Lösung sehr mühevoll ist. Um das Problem möglichst zu vereinfachen, 
approximieren wir die Herrzsche nichtlineare Gleichung für die relative Annäherung durch die 
lineare Beziehung: 
ER (2.4), 


wo x die durch die Einheitsberührungskraft verursachte relative Annäherung ist.!) 


Diese Vereinfachung können wir deshalb machen, weil in den praktischen Fällen der 
Berührung zweier Körper (des Rades und der Schiene) die HErtzschen Voraussetzungen doch 
nicht erfüllt sind (Einfluß der Berührungsflächenbeschaffenheit, der Reibung, der Deformationen 
der Körper auch außerhalb der kleinen Umgebung des Berührungspunktes usw.). Daher stellen 
wir den linearen Koeffizienten x am besten experimentell durch Ausgleichung der gemessenen 
Beziehung zwischen der Annäherung und der Berührungskraft mittels einer Geraden fest. — 
Hertz und RAYLEIGH haben gezeigt (s. [15]), daß die Gl. (2.4) für die statische Lokaldeformation 
(in elastischem Bereich) mit der Gleichung für die Balkenquerschwingung verbunden werden 
kann. Dies wenden wir hier an. 


3. Lösung des Problems 


Unsere Aufgabe wird durch folgende 4 Gleichungen bestimmt: durch die Gleichgewichts- 
bedingung für die Kräfte, die 1. auf die Masse m,, 2. auf die Masse m, wirken, 3. durch die für den 
Balken geltende Querschwingungs-Differentialgleichung und 4. durch die aus der Voraus- 
setzung, daß das Rad stets im Kontakt mit der Schiene (P(t) > 0) im Punkte x = 0 bleibt, 
folgende Gleichung. Also: 


Ba AR. CE Ne. @.1), 

Pa — m, 70) »G [u„(t) u v,(t)] eh [ön(E) i1# d,()] 3 Pd) Nr 8.2), 
EIN Ed Lude,d) FT2E mis d Fk, dd = Pl) ...: (3.3), 
OA I Ba (3.4), 


bei P(t) > 0, samt Anfangsbedingungen für t = 0: 
PX0) —0, 2,0) = 0; 70,0, 2 0 Frl) = rl), 2,0)—=0 8.5) 
und Randbedingungen für |x]— oo: 
ER an DE en Er AN a EN 
Hier bedeutet: 


P, = m,g, P„ = m„g mit der Erdbeschleunigung g, 
P=P,+P,„ den Raddruck; 
v,(i) die Verschiebung der Masse m,, v,(t) = dv,(t)/dt, v,(t) = d2v,(t)/dt2; 
v„({) die Verschiebung der Masse m,, v„(t) = dv„(t)/dt, vd„(t) = d?v,(t)/di2. 


Die Verschiebungen v,(f) und v,„(t) werden von der Gleichgewichtslage im statisch belasteten 
Stande gemessen. Die positiven Richtungen gehen abwärts. 


C die Federkonstante; 
C, den Federdämpfungskoeffizienten [kg em! s]; 

P(t) die zwischen Rad und Schiene im Punkte x = 0 wirkende Einzellast. 
P(t) kann man zerlegen: P(t) = P + R(t). Dabei ist R(f) die dynamische 
Zusatzbelastung. 

E Elastizitätsmodul des Balkens; 
J konstantes Querschnittsträgheitsmoment des Balkens; 
v(x,t) die Balkendurchbiegung an der Stelle x und in der Zeit £ (positiv abwärts). 
v(x, {) wird vom unbelasteten Stand gemessen; 


vlz,t) = Qu(e, Dj0x, vlt) = our, 1)jor”, va td. Que, t)jea®, 
vWz, 1) = Ava, dort, da, t) = Aula, tet, dx. 1) = Ava, N/OR ; 


1) DAVIDENKoV löst in [16] (1949) den Querstoß auf einen einfachen Balken, den er als System mit 
einem Freiheitsgrad vereinfacht, und benutzt dabei auch die lineare Beziehung (2.4). Die Abweichungen seiner 
Lösung von der TIMOSHENKos sind sehr gering. 
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die Balkendurchbiegung, verursacht durch eine in x = 0 bestehende sta- 
tische Last P; v, = vp(0). Von vp(x) wird die durch die Kraft R(t) ver- 
ursachte Zusatzdurchbiegung vz(x, f) gemessen: v(x,t) = Dp(X) + vr, t); 
u die Balkenmasse pro Längeneinheit; 

&;, die Dämpfungskreisfrequenz im Balken; 


Dp(X) 


den Koeffizienten derelastischen WınkLexschen Unterlage [kg/cm?] (s.[9D; 


k 
ö(x) die Diracsche Funktion. 


Vom Gesichtspunkt der Schienenbeanspruchung werden wir uns besonders für die Durch- 
biegung v(x, t), die Kraft Pt), das Biegungsmoment M(z, I) = —EJv(x,t) und die Biege- 
spannung der Randfasern o(x, ) = M(«, !)/W (W — Widerstandsmoment der Schiene) interes- 
sieren. Die für die Tragfähigkeitsbeurteilung entscheidende größte Beanspruchung wird er- 
sichtlich im Lastangriffspunkt entstehen. Infolgedessen reichen wir mit der Berechnung von 
v(0, t) und M(0, t) aus (im Punkte x = 0). Die Frage der Ausbreitung der Durchbiegungswellen 
wird hier nicht untersucht, denn sie würde eine“spezielle Studie erfordern. 


Unser Problem wird mit Hilfe der Integraltransformationsmethode, u. zw. mittels der 
LaPLacE-Carsosschen (auch Wagnerschen) Transformation (Abkürzung L. C. Transf.): 


c-Fioo 


Bild: F(p)=p fee zent, Original: /() = Be j ig ur dp 
a 2ri p 
und der Fourızgschen Transformation (Abkürzung F. Transf.): 
Feg)= [fa)etwdr, Id) = 5 | Br: 


gelöst. Übliche Operatorenbeziehungen werden aus den Transformationentafelwerken [17] und 
[18] ohne Zitierung benutzt. In der Arbeit bezeichnen wir: L.C. Transf. (Bild ): 


Velz,p) =pfvr&,t)erid; Mr@&p)=pf Mr, t)ertdt; V.(p)=p f v„(t) ePt dt; 
(0) 0 0 


V,(p) = p full) er! dt; Rip) =p[ROer'd u.ä. 
0) 
und F. Transf.: j 


Vr(g p) ul Val, p) et: dx; M’%(q, pP) = f Mr(z, p) e?* dr; 
Vr(g) =. vp(x) e"t* de; M%(g) = [ M;&) erde u.ä. 


3.1. Das Rad im Kontakt mit dem Balken 
Wir zerlegen 


P=P+ RO VER ddr (3.7), 
wo vp(x) die Lösung der Gleichung 
EI HK) =ER)P un ar. (3.8) 


unter den Randbedingungen für | 0: vp(x) = vL(x) = va) — pu FR - 
sung der Gl. (3.8) ist bekannt, Na [9]: EL) = BER) = ERLEBEN aratelli AD 


P 
xt) = - e Pla (eosßr x RR ? 
DplX) SEI e (cosßx + sinß|x|); vd, = dp(0) = SE : (8.9), 
wo ß?—= k/(4E J) und FF. Transf. V*%(o) = u F L 
K ) ransf. V’(g) arme... Nach [9] auch 
P 
My.(&) = — e-Pilxl sß: x R :; 
PR) = ze" (os ße —sinfla); M, = M;(0) = TB . (3.10) 
und F. Transf. Mb) = p_ 25 
q* -- 4 D* y 


i 
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"Mit (3.7—8) erhalten wir dann statt des Systems (3.1—4) folgende Gln.: 


BP CO OT GT Ren 3.11), 

— m Ön(t) — € [on — ud] — Cd — I — RD =0 2.2.2... 0819), 
EJ vr!) + wör@, !) +2 u @Drl&,t) + kopla,t) = ER EAST EN) 
vr, HR =)... @.14), 


bei R(t) > — P (Rad stets im Kontakt mit Balken), unter den Anfangsbedingungen für = 0: 
v()=0, NV) =0; - v0, 2,0 0% vu: 0) 0, vr(,0)=0 (3.15) 
und unter den Randbedingungen für |x|— oo: 
REST Ne ee ee (3.16). 
Die Durchbiegung v(z, !) = vp(X) + Dvr(x, t) entspricht dann dem ursprünglichen System 
(3.1—4) samt den Bedingungen (3.5—6). Die größte Beanspruchung herrscht in x = 0: 
v0, =w[l + vr(0, to] ; M(0,9) = M,[1 + Mx(,t)/M,] ; 
0) =o[1+Mr(0,d)/Ml; = Pl4B W). 
Wir führen folgende Bezeichnungen ein: 


=Cm; on — Clm,; = + 0; Do = Ws — OB; 
20 = Gm; 20m = Colmyn 5, VB — pr + Opn; 
ah 7 = H; A=uldEN; 
A=p[(p + op)? + oil; e= Alp + on + opt. 
Dann erhalten wir aus den Gln. (3.11—12) nach der L. C. Transf. mit Rücksicht auf (3.15): 
er 2 - 2 > 5 2 
Ve) = A Fe dt. 


L.C. Transf. der Gl. (3.13) mit (3.15) ergibt 
E J Vi@ p) + up? Vrk&, pP) + 2 u @P Vn&, P) + k Vn(&, p) = ö(x) R(p) 
und F. Transf. dieser Gl. mit (3.16) ist 
E J q* Viq, pP) + up? VXq, P) + 2u @p VRXq, pP) + k Vä@, p) = Rip). 
Aus der letzten Gl.: 


R( 1 
Vk@,p) = A) awed 
Die F. Umkehrtransf. dieser Gl. ist: 


— Ri | 
Vr(, pP) = SaET eretl(cosox +sino|x|), 


aus der man bei x = 0 Folgendes erhält: 
> Br R(p) 
Vx(0, p) = ArET 
Wenn wir Vz(0, p) und V,„(p) in die transformierte (L. C. Transf.) Gl. (3.14) einsetzen, errechnen 
wir aus (3.14): Zug 
— — m o®E 
R(p) = — F(p) —— 
(P) (p) Dip) 


wobei Dp)=xm, AÄ8®EJ+mA+8"EJP+2W,Pp+ &:). Nach der Einsetzung 
von (3.17) in Vz(0, p), Mx(0,p), Vn(p), V,(p) erhalten wir: 


en „N „A | 

Vr00, pP) = — F(p) . ESG en 57 (3.18); 

je — A2@EJ 

M 0, p) = — F(p) a RE ar (3.19), 
= - 808EJ (pe+2@,P+ ©) J 
V.(p) = F(p) — p Br De a (3.20); 


Rn — 80EJ(2 0;,P-+ @2) 
V.(p) = Fip) — ei a NT, >21). 
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F(p) bedeutet in allen Formeln die L. C. Transf. der Unebenheitskurve f(f) (die Existenz 
F(p) vorausgesetzt). Im Falle der Flachstelle gilt nach (2.13): 
® 


Fi) = Bo) ie) = Sa apa ler ae 6 


Der Einfluß der Kontaktdeformation ist in diesem Falle sehr gering und wir setzen des- 
halb in (3.17—21) «= 0 ein. Die Originale (L. C. Umkehrtransf.) der hier abgeleiteten Bilder 
werden mit Hilfe der in Abs. 4 dargelegten Methode gefunden. 

3.2. DasRadverliertKontaktmitdemBalken. Wenn wir das Problem des unrunden 
Rades nach Abs. 3.1. (1. Phase) lösen, kann der Zeitpunkt t, eintreten, in dem P(t) = P + R(t,) 
— 0; dann verliert das Rad den Kontakt mit dem Balken. Darauf folgt die Eigenschwingung des 
Balkens und des Systems von 2 Massenpunkten (2. Phase). Die 2. Phase dauert bis zum Zeit- 
punkt t,, in dem das Rad den Balken wieder berührt und stößt. Dann schwingt der Balken samt 
dem System (3. Phase). Die ganze Erscheinung kann sich mehrmals wiederholen, bis das Rad auf 
dem Balken wieder stabil ruht. Deshalb ist dieser Fall komplizierter als der vorige (Abs. 3.1). 
Um ihn zu vereinfachen, berücksichtigen wir, daß die Zeitdauer Z, in praktisch berechneten Fällen 
im Verhältnis zu T sehr gering ist (Ordnung t, =» T/100) und die Kraft P(t) in dieser Zeitdauer 
vom konstanten Wert P sehr schnell auf Null absinkt. Die weitere Phase wird durch die Ent- 
lastungszeitdauer nicht wesentlich beeinflußt. Wir wollen daher die ganze 1. Phase auslassen und 
voraussetzen, daß im Zeitpunkt t = 0 die plötzliche Entlastung (P(f) = 0) eintritt. 


3.2.1. Eigenschwingung — 2. Phase. Die Schwingung des Systems und Balkens wird wieder 
durch die Gln. (3.1—3) bei P(t) = 0 samt Bedingungen (3.5—6) gegeben. Die Gl. (3.4) gilt 
hier nicht. Die Lösung der Gln. (3.1—2), in denen P(f) = 0, mit (3.5) in Bildern ist ähnlich wie 
in Abs. 3.1.: 


= Erde. zu 2 
V„(p) = P ; V.(p) = P re 
Das Original dieser Bilder kann man direkt hinschreiben: 
a 2 A 
DD 7, I. $) F 20, HF 7 (3.23), 
PD 2 
D(= 29 2.4 EROT., Fr ee (3.24). 
n 
Auch die Geschwindigkeiten der Punkte m,, m, sind wichtig: 
2 2 R 
d,(t) = PB ls +2 HEND - - --.... BB), 
: Pr = 
DO = 59 [2 dr Ya) FF aan a * 
n 
wobei die Bezeichnung 
1 -Öönt .: 2 
lt) = — ER sinwgt; 
Dogg z 
It) = a 1— e "3 (cos ost +2 sin w‘ 
p} DRs N dos 4 O8 sın ot S 
1 2 ’2 ® 
een Er \ | 
a er > EUR, 
got) == P} : n 211 ren e”®B" cos et - _e@pt @B 3 005 — WE -Wpt ; ’ 
@s| 2 Or ©0s os Bu ls 2, Sau 1 


benutzt wird. 
L. C. Transf. und F. Transtf. der Gl. (3.3), in der P(t) = 0, mit Bedingungen (3.5—6) ergibt: 


V*(q,p) = P ap+2uo,p 1 h 
EJ .At+4ß 'EJ +4 Ve) + VRo(g P) 


ON IREEEN "USER EEE Lues 
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fi ei 2 7 f En = 
wo Vk,(, pP) = — EI er « F. Umkehrtransf. ist: V (x, p) = vp(x) + Vro(&,p) und p 


=> = 3 
bi 2=0: V(, pP) =W+ Vo, pP) = a, 2 


Bi Ip + o)° + a) Re Un 


transf. dieser Gl. ist: 


BIO: DER, DD AR (3:27), 
wobei 

een 

VRo(0, 1) = Dy (=) ne SEP Idıza 3) de An (3.28); 


weil gilt: 3 


See ü 
a Eh 


P 
[(p + a)? + pep+ir 


für Rev> — Es ‚ md z=0f; !=or; 4d= a. Zerle)* 
[J,(&) — Besseusche Funktion erster Art, Ordnung v; /'(x) — Gammafunktion.] 
Das Biegemoment errechnen wir ähnlich: 
NORDENS AM ElOEL) ae ee), 
wobei 


1/2 9-1/4 T (3/4 a N 
0 m, (2) u ed ee (3.30). 
0 


Die 2. Phase endet im Zeitpunkt i,, wo das Rad wieder den Balken berührt. Den Zeit- 
punkt 2, erhalten wir aus der Gl. (3.4), bzw. (3.14), die hier wieder gilt: 


Dr(0, 13) = URo(O, 12) = # P + vultı) — I); (Be) Der aa 


Wir lösen die Gl. (3.31) am besten graphisch. Der Schnittpunkt der Kurven an der linken und 
rechten Seite (3.31) gibt für unseren Zweck den Zeitpunkt t, genau genug an. 

Die Eigenschwingungsfunktionen (3.27—30) sind numerisch errechnet (0, = 0 und @/oy 
— 0,15) und in Bild 3 gezeigt. Die Eigenschwingung des unendlichen Balkens auf elastischer Un- 
terlage ist sehr interessant: die Energiezerstreuung tritt auch im Falle fehlender Dämpfung ein.?) 


-08 


08 :3.00:4. &r=015 


Bild 3. Eigenschwingungsfunktionen 
(0,2) M(0,t) 


:3.0,=0; 4. 0,j0 = 0,15) 


:1.9,=0; 2, 2/0, = 0,15, Biegemoment 


r D) 
(Durchbiegung M, 


2) Auch Esmeiser [19] (1948) löst die Schwingungsaufgaben des unendlichen Balkens auf elastischer 
Unterlage, gibt aber nicht die Formeln für die Eigenschwingung an. (Von dieser Arbeit erhielt der Verfasser 
erst nach der Beendigung des Manuskripts Kenntnis.) 
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3.2.2. Stoß — 3. Phase. Die 3. Phase beginnt im Zeitpunkt t,, wo die Zeit in dieser Phase ihren 
Anfang nimmt. Die Anfangsbedingungen für diese Phase sind folgendermaßen gegeben (für 


_ ” v,(0) = Dy2 > d,(0) = (735 v„(0) = Da» d„(0) a En . (3.32); 


v0) =), RN) =ER) ..-- ernennen (3.33). 
Dabei stellen 
Dr = lb) » Ga = D,(t,) 5 Ung = Vn(t,) , CH = D„(t) ; 
v(X) = vp(X) + Vrol& bo) » C2(2) = Urol&; 1;) a 


die Werte am Ende der 2. Phase dar. Die Randbedingungen werden wieder durch (3.6) bestimmt. 
Die GlIn. (3.14) und mit Hilfe von (3.7—8) auch die GlIn. (3.11—14) lösen unsere Aufgabe unter 
den Anfangsbedingungen (3.32—33). Aus (3.11—12) und (3.32) erhalten wir dann in Bildern: 


20H P+ © 


m, A 


+ B(p), 


te 7; Vi =— Rp) 
wo r 
A(p) = Up? + 2 0, P + ©) [(p? + 2 @yn P) Una + P ine — 2 @bn P Vra] 
+ (2 @nPp + or) [(P? +2 @rP) Ura + P ra — 2 @yr P Una} ; 
Bio) = Up? + 20m p + @B) IB? + 20m P) ra + Pa —2 0 Pina] 
+ 2 0, pP + @) [(P? + 2 @pn P) Una + P in — 2 @®pn P lra]} - 
Die Originale zu diesen Gleichungen sind: 
AM = n2 9a) + (2 ©B Ung + Cn2) Is(t) + (OF Ung + On Urn + 2 Opr Enz + 2 @pn Cr2) galt) 
+ [®r ng + @n Cr2 + 2 (pn ®r — Opr On) Ung + 2 (Opr © — pn ©) Dre] 9ıld) ; 
Bl) = 1,5940) + (2 @> dr + Cr2) Iz(t) + (On Org + ©F Una + 2 @pn Cra + 2 @pr Enz) galt) 


+ [@ 2 + @r na + 2 (@,.@ — OO) Ua 2 (Opn © — Opr 05) Oma] 9 
wo 


-Opt DBi% ‚ 
gl) =Ee ? (cos gg E— —- SIN @gs ) : 
os ‚ 


. Für den Balken gilt die Gl. (3.13) mit den Randbedingungen (3.16). Die Anfangsbedingungen, 
mit Hilfe von (3.33) dargestellt, sind: 


vr(X, 0) = Vr5lt) , Dr(%, O) = cpul£) 2 Mor (3.34) 


wo Urg(X) = Vrol&; La), Cra(X) = Uro(x, t), was aus den Endwerten der 2. Phase folgt. Dann wird 
die L. C. Transf. der Gl. (3.13) samt (3.34) folgendermaßen aussehen: 


EJV!(«,p) + up? Vale, pP) + 2u@p Vr(a, P) + kVplc, P)— (up? +2 u 0, P) Urz(X) — UP Cre(X) 


| — ö(2) Rip) ; 
die F. Transf. dieser Gl. wird mit Hilfe der Beziehungen 
Vre(g) = Vhoq L)=— Si . 11er cost eine ll 
2 EJ gi+4ß "RER T In, EN 
* = 271 N ! 2 Ay) f 2 P 
VRo( pP) = p J Vro(g tb) err'dt; = Pla +4P) — w,; P=E Jju 
und 
ER T a l > 
Che(q) = VR 2) = — a EL NWE 
Ra(q) Ro(9; l,) E 7 gi 2 4 pi (e E= rg e b sın 0% t, 


die Form bekommen: 
EJg' Vin pP) + mp? VRQ, pP) + 2u 0.p VA, p) + kVklq, pP) — (up? + 2u@,p) Vho(gst,) 


— up VRo(b) = Rp). 


en 


- 
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Daraus folgt: 


R(p) 1 2+2 : 
Velo, D) — np? + 2u0,p up = 
r(q P) E 7 gt 2 4 ot EJ (q! ns 4 0%) VRo(q t,) + E 7 (g* _ 4 1) Vro( t,) £ 


Die letzten zwei Glieder dieser Gl. sind aber die einfachen Bilder der bekannten Funktionen: 


Rp) 1 
Vr(@, pP) = ET q + 4 ot ; + VRoa(g; pP) — VR0(g P) 


wo 
Urox(& E) = Urn & +1) 
und 
z 21 2uw P 1 
verap)o-ve Ep rawoyp IR = plz Ur 
02(9: P) SE sagen wwer, EJ a un sin o, t, 


Kur P 1 A 
BETTEN "EJ m+4p en sin o,b%,. 


Bei x = 0 (v2o(0, t) ist durch (3.28) gegeben): 


el Devon) 


soEJ 
Wir setzen Vz(0, P); V„(p) in die transformierte (L. C. Transf.) Gl. (3.14) ein und erhalten daraus: 
— = mASeERJ 
a ee 3. 


wo = 


cp) = Afp) — Vro(0; pP) + Van, Rp; For fhWerd; KO=il+Dd, 


und weiter 


Vn(0, p) = Ep) 7 : ; De een (3.36); 
1400, p) = Co) ET + Mal, Mn) 83: 
Ver) = N a... (3.38); 
ec nn u EN ....@.39), 


wo — ähnlich wie früher Mar95(&; i) = Mr0(&; t; + ) und Mzo(0, £) durch die Gl. (3.30) gegeben 
wird. 

Die Originale findet man nach Abs. 4. Der Koeffizient x hat beim Stoß einen größeren Ein- 
fluß besonders auf den Kraftverlauf. Darum vernachlässigen wir ihn in diesem Falle nicht. 


4. Numerische Berechnung der Laplace-Carsonschen Umkehrtransformation 


Es hat sich als vorteilhaft erwiesen, unsere sehr komplizierten Bilder in solche Beziehungen 
umzuformen, welche die numerisch leichter lösbaren VOLTERRASchen Integralgleichungen dar- 
stellen. Der Verfasser hat die Methode von PETERKA [20] (1956) auch auf die VOLTERRAschen 
Gleichungen erster Art erweitert und die Lösung mit Hilfe der Smmpsoxschen Regel präzisiert.?) 

Der Grundgedanke der numerischen Methode besteht darin, das Bild Y(p) der gesuchten 
Funktion y(f) auf folgende Form zu bringen: 


— ,_ U) | 
ee ae 
che 

bzw. 2 
= U(p) ’ 
Y BE LIE Mr ee ee (AD): 
zoo nn 


3) Die eigentliche numerische Lösung der VorLrerraschen Integralgleichungen (zweiter Art) ist schon 
lange bekannt, z. B. MikeLAnze [21] (1935). Ich weise auf die Fehlerabschätzung in [21] hin. 
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wo Up), Z(p) die Bilder sind, deren Originale u(), z(f) bekannt sind oder leicht berechnet werden 
können. Das Original der Gl. (4.1) bzw. (4.2) ist: Bu - 


fx C-Dyda=ul) ......h.rn: (43), 3 


bzw. 


y() + f ze -NDyd)d=uh) ui onen. (4.4). 


Das sind Votrerrasche Integralgleichungen erster bzw. zweiter Art. Das Faltungsintegral in 
den Gln. (4.3—4) kann man numerisch nach der Smmpsoxschen Regel (oder nach jeder anderen 
Methode der numerischen Quadratur mit äquidistanter Teilung) unter Zerlegung des Intervalls 
0 bis tin j Teile der Größe h integrieren: 


h . } 
Für gerades 7: 3 3 +4Ay3-1 +2%23+'+24 32 +4, 14 +Y2)- 


Für ungerades j [die gerade Teilzahl besteht im Intervall O bis (/ — 1) h; für die Berechnung 
des Integrals im letzten Intervall ersetzt man annähernd die Kurve z (’—r) y(r) wieder durch eine 
durch die Punkte j — 2, j— 1, j bestimmte quadratische Funktion (Parabel)]: 


h 15 B) 
(mat ins t2mas th +2 n + g Han +3natqun), 


wobei y;, z;, u; die Werte der Funktionen y(j h), z(j h), u(j h) darstellen. Dann kann man den un- 
bekannten Wert y; mit Hilfe der früher errechneten y,, Yı, - - - , 4;_] bestimmen. Das Verfahren 
kann tabellarisiert werden. Die Tabelle für die numerische Berechnung der VoLTERRAschen Inte- 
gralgleichung erster Art (4.3) sieht bei = 3u;/(hz,); 2; = z;/z, (in unseren Fällen z, — 0) 
folgendermaßen aus: 


7, = lim y(f) = lim Y(p) ) 
0 Pp>X 


18 2 
Y=7W%— 2) 


Ir/& pe De 
y=7 BT re 
Ter® Ei is. wr 
Yg— A (u — yYy2—4yı 3 — 2002) 3 Wu, Ver (4.5) 
ee 3 = == re 
Yı= (us Ina 2a na) 
1 


I y % Yo Ze 4,3 —2, 3 —4y 233 —2y 23) 


Yy=Wı 


h h 
Yı == U, os Sg Yo 2, 1 nr 6 a) 


h 
»=u7hr+4yz) 
7 h 15 } 
a ee 3 (402 E. zYı% +31 ) us, al (E 


h 
Ku Yattyza +23 +Ayz) 


Fr h 15 
„eu 3 (u 3 +4y24 +2 Ya + eu Yy2y,+3y a) 


tzung der Funktionen z (l— 7) und y(r) durch 


(Der Wert y, wird mit Hilfe der annähernden Erse 
die Geraden errechnet.) 


ür den Fallx — O0 zweckmäßig. Für x + Oist dagegen 
er numerischen Methode: 1. Verschiedene Möglich- 


#1. DasRad im Kontakt mit dem Balken. Die Bilder der gesuchten Funktionen sind : 

in 8.17—21) gegeben, wo x = 0. Wir bringen sie mittels Division des Zählers und Nenners durch 
_ den Ausdruck m, A 80?EJ auf die Form (4.1). Dabei benützen wir mit Vorteil dimensionslose 
i ‚Parameter und Veränderliche £ = wt, &= wor. Wir bestimmen auf diese Weise die Funktionen En 
2), u(&), diein Taf. 2 angegeben sind. l 


Tafel 2. Funktionen z(&) und u(z) 


= wt 3 
!=or 2@) 


A ETF NE u END 


eo SE Er a a + | 
2 
1 + B (ev — @B a) ( a BEN cos = ; — a gBzı®, sin es Pi er 
2 nn nn nen > 
| u(#) 
R(&)  - 
4 PB San sin = (4.8) 
a0.) ee En Fra On [ra 000 £as 


[O7 ©, 


ae kin gar) (4.9) 
0 4 


E 40% a/o\tt /oo\ 1 2ATT)[E m 
=) (&) Tee 24,8) dE 
9) & 


2 \ 2 um 250 Ar = en 
ge & — | cos& + er ® sin 2 ! 2d,,,(£) cos & dE 


D 


25 . - 2 & = - u Shot 
+ 2 ) sind — 03 cos 2 1 2d,,,(&) sin & |, | (4.10) 


E E > 
wo 2d,(&) = e’ ER ‚ze, J, (0, &/®). Die Funktion y(&) wird dann direkt aus der Gl. [z 3) 
0 


x y(&) dE = u(&) errechnet; für £ > oT findet man die zu bestimmende Funktion mittels der Be- 
ziehung y(£) — y(£— oT). Die Funktionen u(t) werden ebenfalls numerisch (nach Simpson) 
G+2)h 
z h \ h 
berechnet: u;,, = j o(E) de — u; + 3 (9% + 49,1 + Piz) > dabei a ie 
0 


4.2. Das Rad verliert Kontakt mit dem Balken. Die Lösung der 2. Phase wird in 
(3.23—30) gegeben. Für die 3. Phase gelten die Bilder (3.35—39), wo x =# 0. Division des Zählers 
und Nenners der Bilder durch den Ausdruck x m, 4 8 0? E J führt auf die Form (4.2). Im Falle 


des Stoßes benützen wir die dimensionslosen Veränderlichen € = wpt, &= or. Die Funktionen 
2() und u() sind in Taf. 3 angeführt. 


4) Die VoOLTERRAsche Integralgleichung erster und zweiter Art wurde in [1] auch mit Hilfe der Reihen- 
entwicklung gelöst. Diese Methode ist vorteilhaft für kleine Argumente, für die größeren dagegen ist die Reihen- 
konvergenz zu langsam. 
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Tafel3. Funktionen z(z) und u(z) (für 3. Phase) 


Bol 2(®) - | 
E= oT n a — 
v, (o\MT(5/4),, 
|) Ta Re 
1 [Dy; ER N Se. ©2/®,, in in 98 z TEN 
e: u Mn 05° Er 3 08 Or @i8 22 ” 
, \ ee 2 0) @br @B -og2/o wS- 3] 
WFT IE - opE/® 8 1 7 —e -sin — (4.11) 
ano Rt , BE/o, 0. ea 2) og Pop o, $ 
mn m ne 
z =. _ 
y@) ua) = [y@— 3) e@) ds | 
[) 
= 0(&) abe (4.12) 
0: 0 OB... Jaja [oeytle BaFin 
en x) Y(%) = =F v% 3 PX) — () T (3/2) 114(%) Yo 
vr(0, x) _ vRos(O,t) vro(O,t), wr(0,0) _ wRoz(0,0) _ vRolO, to) - ' (4.13) 
% = y(%) Hr © % > Vo © U N 
— Ö(R) für WELCH 2. 60) 
M 0, — Az . e eFplzi Self . 
a z er nl Eee! I" a NER 
“Pu, 
g ya nn (3/4) __ Mro(0, t) _ gen ni 
9%) -(®) 5 T (1/2) f 1d_ 1.(E) dE M ’ = 0 ’ Yo 03 | 
MO) _ zen, Mrd) MrolOi)  Mr(0,0) _ Mros0,0) _ Mr) | (419) 
U Ag, He SUR M, Mr, | 
dc@) _ 1 dc) 1 dA or ETC), ar 13 
Ne Ta aaa) 
a 212 7 (5/4), RO Em. sin - Ran. 
e) T (3/2) ao o, (en Bu 


z z > = .- 
Die Funktion y(®) wird numerisch aus der Gl. y@) + fz (€ — y(d) d£ = u(t) errechnet. Für 
0 


die Funktion u(«) gilt in diesem Falle das Faltungsintegral: u) = j y (© — &) olE) dE, das nach 
SIMPSON folgenderweise tabuliert werden kann (9; = go(j h), y; = y ( h), bi a, = 0): 


w=0 

h | 

3 (Yo Pı 5 Yı9ı 
h 

op + Ay) 


h 15 1 
-3|%Pr ua 3 Ya Pı 


a 


ja‘ 


(4.15). 
h € 
RT (Yopı +4 Yı9 +29 +4 Y;9ı) 


>. Vergleich der Theorie mit dem Versuch 
Um die Richtigkeit. der Theorie zu zeigen, haben 
durchgeführt, die in [1] näher beschrieben ist. "Das Modell 


wir eine Versuchsmessung am Modell 
lischen Ähnlichkeit. Da die theoretische Lösung v 


erfüllte nicht alle Gesetze der physika- 
orhanden war, war das auch nicht Rue Es 


E 
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wurden 4 charakteristische Fälle herausgegriffen, nach der Theorie eingehend berechnet (Zeit- 
verlauf der Kraft, Spannung und Durchbiegung) und mit dem Versuch verglichen. Dabei hat 
sich gezeigt, daß die Theorie mit der Versuchsmessung sehr gut übereinstimmt. — Wir haben 


_ vorläufig keine eigenen Ergebnisse von Messungen an den Schienen in der Naturgröße und des- 


halb haben wir außer der Messung am Modell auch die publizierten Ergebnisse der amerikani- 
schen Messungen [11] zur Überprüfung der Theorie benutzt. In [1] wurden 3 charakteristische 
Fälle des theoretischen Schienenspannungsverlaufes errechnet, deren Experimentalverläufe in 
[11] publiziert worden waren (Fahrt Nr. 238, 252 und Fahrt bei der kritischen Geschwindigkeit 
20 Meilen/St.). Auch hier wurde eine gute Übereinstimmung festgestellt. 


ai 


N 


Bin 
EBEHEN 
ERBEERENISG! 


Theorie S=-SBI9Klem®———_—---—- Experiment 


Bild 4. Spannungsverlauf bei kritischer Geschwindigkeit (die experimentelle Kurve 

ist aus [11], $. 440, Registr. 20 mph, 100 RA-A Rail, entnommen, u. zw. nach 

Abrechnung der Spannung, welche dasselbe Rad ohne Flachstelle verursachen 
könnte) 


Bild 4 bringt als Beispiel den berechneten und experimentalen Spannungszeitverlauf für 
die Fahrt bei kritischer Geschwindigkeit (die experimentelle Kurve aus [11], S. 440). Bei der 
Berechnung wurde die Eigenfrequenz des Oberbaus n, = w,/(2 x) nach [11] und besonders nach 
den Messungen von BETZHOLD [13] als der Wert n, = 60--70 Hz eingeführt. Zu diesem Wert 
gehört in der Formel & = (2.n,)? = k/u eine wesentlich größere Masse u, als die Summe der 
Massen von Schiene, Kleineisen und Schwellen beträgt. Das bedeutet, daß hier teilweise auch 
die Unterlagmasse mitschwingt. 


Die zur Berechnung benutzten Werte wurden für unser Beispiel auf Bild 4 zum Teil aus dem 
amerikanischen Bericht [11] entnommen (diese Werte sind im Folgenden ohne Bezeichnung an- 
geführt). Da aber dieser Bericht leider nicht alle zur Berechnung notwendigen Angaben zur Ver- 
fügung stellt, wurden weitere aus dem Spannungsverlauf errechnet (im Folgenden werden diese 
Werte mit * bezeichnet), und die Angaben über den Wagen wurden aus [23] von den für das 
Fahrgestell ,„AAR Truck“ für Güterwagen festgestellten Werten entnommen (und im Folgenden 
mit ** bezeichnet). Die zur Berechnung gebrauchten Zahlenwerte sind also: unrundes Rad 4!/, Inch, 


ad = 0,47cm; mW = 82,4 st; ' al — lo 

We 125729}; ww =; Mir — 20 Meilen/St. 

P =11480kg; Schiene 100 RA—A: — 32,2 km/h = 894 cm/s; 
Pr = 085 kg**; E = 2,1 10° kgjcm?; ’“ —= 1,423 : 1025; 
Pz0795 ke **, I -=i2.100 cmf; Dee 

C = 4460 kg/em**; W = 254,4 cm?; g/(v, ©) = 0,0978; 

12 ALM; ra = 069 kölcher, (g P, ®)/(v, @ Pn ®os) 
RR Ba DA; — 8,26; 

Be, ; ee NE le Fehr | osl® = 0,1867. 


Das Biegemoment wurde aus diesen Werten nach den Gln. (4.7), (4.10), (4.5), x =0; 
h=2n/10;j=0,1,2,...,20 errechnet. — In diesem Falle verliert zwar das Rad den Kontakt 
mit der Schiene, aber die Stoßgeschwindigkeit ist gering, so daß wir angenähert auch mit negativen 
Kräften zwischen Rad und Schiene rechnen können. Dies zeigte sich auch bei Modellmessungen 
in [1]. — Aus dem Bild 4 ist ersichtlich, daß das maximale positive Verhältnis Mx(0, DM, = 2,75 
beträgt; die maximale Biegungsspannung in den unteren Fasern der Schiene ist also: 


max o(0, {) = o, [1 + max Mx(0, {)/M,] = 589 (1 + 2,75) = 2 210 kg/em?. 


ih 


Der Einfluß der Flachstelle auf die Schiene wird für 
maßen idealisiert: die Schiene als ein unendlicher, auf einer elastis 
und der Eisenbahnwagen als ein System mit 2 Freiheitsgraden. Das ur a 
Stelle (ohne Reibung) drehen. Die Schwingung des Balkens wird durch lie . 
zwischen Radschwerpunkt und Balkenmittellinie erzwungen. Die nic line 
schen der Berührungsverformung und der Berührungskraft während des 5, 
durch eine lineare Beziehung ersetzt. Es wird ein System von Differentialg] 8 
das durch die Integraltransformationsmethode gelöst wird. Die LAPLAcE-CARsoN 
werden auf Formen gebracht, deren Originale die Vorrerraschen Integralgleichung stellen; 
diese werden dann numerisch gelöst. Die Lösung ermöglicht es, den Zeitverlauf der Kraft, S ei 
nung, Durchbiegung, die Zeitdauer des Stoßes und andere Werte zu errechnen. Es wird auch n 
numerisches Beispiel berechnet und das Ergebnis mit einer Versuchsmessung verglichen r 
Lösung gilt für die Unregelmäßigkeiten sowohl an den Rädern (z. B. Flachstellen), als auch auf 
der Fahrbahn (z. B. vereinzelte Schienenunebenheiten). Sie ‘kann auf andere Konstruktionen 
sowie auch im Lokomotiv- und Wagenbau übertragen werden. g 


Aus der theoretischen Untersuchung sowie aus der Versuchsforschung [1], [7], [11], diedie 
vorgelegte Theorie bestätigt, folgt, daß die Biegungsspannung in unteren Fasern der Schienen 
sogar um ein Vielfaches größer sein kann als die durch die statische Belastung verursachte Span- = 
nung. Auf Bild 4 beträgt diese Vergrößerung das 3,75fache. In der Praxis schwankt sie zwischen 
1,0 bis 4,5 (s. [11]). Die Schienenspannung ist durch 5 wesentliche dimensionslose Parameter 


bestimmt: 


a {7} g gen Per 0 ae 


; "gs ER! Tr 
Dy [O7 KON) v0 Pn @os {0} 


wie es aus den Gin. (4.7), (4.10) folgt. 


Um den Einfluß von einzelnen Faktoren zu zeigen, haben wir in [1] die Spannungsverläufe 
für einige wesentliche dimensionslose Parameter berechnet. Die Berechnung hat gezeigt, daß die 
maximale Biegespannung mit zunehmender Geschwindigkeit bis zu einer bestimmten kritischen 
Geschwindigkeit steigt und danach ein wenig sinkt. Dies stimmt völlig mit den Versuchsmes- 
sungen sowohl am Modell, als auch in Naturgröße überein, s. [11]. Diese kritische Geschwindig- 
keit beträgt nach den Versuchen auf Schienen [7], [11] etwa 30 km/h. 


Zum Schluß möchte ich den Herren Prof. Ing. Dr. KoLoUSER und Ing. Dr. BABUSKA für 
wertvolle Ratschläge zu meiner Arbeit und den Herren Prof. HEINRICH und Prof. MARGUERRE 
für die Bemerkungen, nach denen ich mich bei der endgültigen Formulierung des Manuskripts 
gerichtet habe, meinen besten Dank zum Ausdruck bringen. 


Literatur 


[1] L. Frysa, Vliv plocheho kola na nekoneönd dlouhy nosnik na pruänem podklad®. Prag 1958. (Diss., 
TH (CVUT), F. f. Ingenieur-Bauwesen, Prag — unveröffentlichte Lithographie). 
[2] H. SarLLer, Zbl. Bauverw. 41 (1921), S. 115—118, 126—129. 
[3] S. Timos#enko, Proc. Int. Congr. Appl. Mech. Zürich 1926, S. 407—418. 
[4] A. Anter, Organ Fortschr. Eisenbahnwes. 90 (1935), S. 173—175. 
[5] Techniteskij spravo&nik Zeleznodoroänika. Bd. 5. Moskau 1951, S. 234—237. 
[6] C. Porr, Arch. Eisenbahntechn. 1 (1952), S. 11—26. 
[7] H. Russ, Arch. Eisenbahntechn. 1 (1952), $. 28—32. 
[8] P. I. Gor’kov, Trudy mosk. inst. in2. Zeleznodoro2. transp. Bd. 80/1. Moskau 1955, S. 189—212. 
[9] L. Frysa, Aplikace mat., Prag. 2 (1957), S. 105—132. 
[10] H. Sarer, Organ Fortschr. Eisenbahnwes. 88 (1933), S. 183—188. 
[11] Bull. Amer. Rly Engng. Ass. 53 (1951), 8. 423—448. 
[12] H. Rusın, Arch. Eisenbahntechn. 2 (1953), S. 51—54. 
113] Ch. Berznouo, Glasers Ann. 81 (1957), 8. 76—82, 108—115, 137—145. 
[14] Handbuch der Physik. Bd. VI. Berlin 1928, Kap. 3., Abs. 79, 8. 271— 276. 
[15] S. Trmoscnenko, Z. Math. Phys. 62 (1913), S. 198— 209. 
[16] N.N. Davıpenkov, Sbornik trudov inst. stroitel. mech. Akad. nauk USSR, Nr. 11 (1949), S. 73—82. 
[17] V. A. Dırkis u. P. I. Kuzw&oov, Spravoönik po operacionnomu iscisleniju. Moskau 1951. 
[18] G. A. Campgerr, u. R. M. Foster, Fourier Integrals for Practical Applications. New York 1948 
[19] W. L. Esmeiser, Appl. Sci. Res. A1 (1948), S. 151—168. 
[20] Souhrn praci o automatisaci 1956. Prag 1957, S. 159—182. 
[21] 8. E.: Mikntadzs, Izv. Akad. nauk SSSR, OMEN. VII (1935), S. 255 —300. 
[22] Tables of Bessel Functions of Fraetional Order. Bd.1. New York 1948. 
[23] Tests of Truck for High-Speed Freight Service. Office of Mechanical Engineer, April 1940. 


Manuskripteingang 26. 11. 1958 


Anschrift: Dipl.-Ing. LAnısLav FRrYBA, Prag 12, Stalinova 67, ÖSR 


‘ H. NEUBER, Synthese der linear-elastischen Schwingungssysteme — Kleine Mitteilungen 185 


ZAMM 40 (1960) Heft 4, Seite 185 


u u a 


Synthese der linear-elastischen Schwingungssysteme”) 
| Von H. NEUBER = 


4 Durch Einführung der inneren Kräfte anstelle der Verschiebungen als Unbekannte ergibt 
sich für die Bestimmung der Eigenfrequenzen eine Determinante von wesentlich einfacherem Auf- 
bau. Die Federkonstanten treten nur noch in den Diagonalfeldern und dort nur im Nenner des 
Frequenzquadrates auf, während alle übrigen Glieder lediglich die Kehrwerte der Punktmassen 
enthalten. In Verbindung mit Regeln der Algebra läßt sich zeigen, daß sämtliche Koeffizienten 
| der Frequenzgleichung auf die Eigenfrequenzen von Ersatzsystemen zurückgeführt werden kön- 
nen, welche entweder durch Erstarrung oder Herauslösen einzelner Federn gewonnen werden 
können. Hieraus ergeben sich neue Formeln für die Berechnung von Frequenzschranken, welche 
auch für Kontinua gelten. Als Beispiel wird gezeigt, wie sich für Dreh-, Longitudinal- und Biege- 
schwingungen von Stäben bestimmte Integrale angeben lassen, welche neue Frequenzschranken 
liefern. Die ausführliche Behandlung der Theorie ist im Aufsatz ‚Über allgemeine Eigenschaften 
der Schwingungszahlen linear-elastischer Systeme‘, Ing. Arch. 28 (1959) S. 229—241 (Festschrift 
f zum 70. Geburtstag von RICHARD GRAMMEL) dargelegt. 


Anschrift: Prof. Dr.. H. NEUBER, Starnberg b. München, Max-Emanuel-Str. 16 


*) Nachtrag zum Tagungsheft ZAMM 39 (1959), H.9—11: Auszug aus dem auf der GAMM-Tagung 1959 
(Hannover) gehaltenen Vortrag des Verfassers. 
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Sequenzanalyse für zwei Alternativfolgen 


1. Einleitung 

Bei einem Experiment seien zwei stochastisch un- 
abhängige alternative Zufallsfolgen I(N,) und II(N,) 
von den Umfängen N, und N, entstanden, deren 
Merkmale durch + bzw. — gekennzeichnet werden 


sollen: 
I: +— ++ -—...N, Elemente (1) 
I + ++ =...X, Elemente > 


p, und 9, seien die unbekannten Wahrscheinlichkeiten 
für ’+’ in den zu I bzw. II gehörenden unendlichen 
Grundgesamtheiten, auf die zurückgeschlossen wer- 


den soll. , =1—p, und, =1-—.p, sind die ent- 


sprechenden Gegenwahrscheinlichkeiten für ’—'. 


(Falls die Grundgesamtheiten endlich sind, dann sind 
?1 Pa 9 und 9, die entsprechenden relativen Häufig- 
keiten.) Welches p,— p, läßt sich in dem durch 
N = Min (N,, N,) charakterisierten Versuchsstadium 
mit dem Aussagerisiko « sichern ? Die Aussagesicher- 
heit 1—« gibt die Wahrscheinlichkeit an, mit der 
das Verfahren zu einer richtigen Entscheidung führt, 


Bei geeigneten Idealisierungen tritt diese Frage 
z.B. in der Medizin auf bei der Entscheidung der 
- Frage, welche von zwei Behandlungsmethoden für die 
gleiche Krankheit die bessere ist, in der Technik bei 
der Bestimmung der günstigeren von zwei Fertigungs- 
arten usw. Die klassischen Methoden zur Absiche- 
rung des Unterschiedes der Mittelwerte zweier Bino- 
minalverteilungen erfordern i.a. numerische Rech- 
nung, die der Praktiker oft scheut. 


Der hier im Anschluß an die Methoden der Sequenz- 
analyse nach A. Warp [4] erstellte vereinheitlichte — 
d.h. für alle « und alle p, — p, gültige — Folgetest- 
plan!) erlaubt in jedem Versuchsstadium ohne vor- 
herige Schätzung von ?,, p,, ohne Fallunterscheidung 
und ohne asymptotische Betrachtungen die graphi- 
sche Ermittlung der möglichen Aussagen. Man er- 
kennt leicht, wann es sinnvoll ist, einen Versuch abzu- 


ı) Für Anwendungen vgl. [1] 


brechen, und ob sich etwa während der Durchführung 
unerwünschte systematische Änderungen der Ver- 
suchsbedingungen eingeschlichen haben. 


2. Durchführung 
In den beiden Folgen (1) treten die übereinander- 


stehenden Zeichen in den Kombinationen (7% (E) A 


() (&) mit den Wahrscheinlichkeiten 9, 95, 91%» 


9ı Ps, Pı 9 auf. Die ersten zwei Kombinationen tragen 
zu Differenzaussagen nichts bei, wir sammeln sie nur 
in einer Strichliste auf. Im Versuchsstadium N mögen 


P(N) (#)-Kombinationen und M(N) (Z)-Kombina- 
tionen vorliegen. Die im Test verarbeiteten verblei- 


benden Kombinationen treten nach der Teilungs- 
regel mit den Wahrscheinlichkeiten 


Re Qı Pa 
Een %t UP 
auf. Die erste Relation zwischen 9, und 9, mit p,_ 
als Parameter liefert Tabelle l und Bild l. Ne) 


Bei bekanntem p, legt p,_, eine bestimmte Differenz 


bzw. P (+) —al 2, (2) 


Pa — p, fest, insbesondere kennzeichnen P(-) >0,5 


bzw. = 0,5 bzw < 0,5 der Reihe nach 9, > p, bzw. 
p, = pı bzw. p < pi. 

Ausgehend vom Ursprung eines xy-Koordinaten- 
systems zeichnen wir fortlaufend mit der Entstehung 
der Folgen I und Il einen Treppenzug, und zwar gehen 
wir, jeweils im Endpunkt des bereits gezeichneten Po- 
lygonzuges beginnend, beim Anfall einer Kombina- 
tion (+ um eine Einheit noch oben, bei En) um eine 
Einheit nach rechts (s. Bild 3). Insbesondere haben 
gı = bzw. 9, = 0 wegen p _, = 0 eine Wanderung 

+ 
längs der x-Achse zur Folge, analog bedeutet q, =1 
bzw. ,=]1 mit pP.) = 1 Wanderung längs der 
. 
y-Achse (vorausgesetzt, daß nicht nur Kombina- 
tionen (7) bzw. (=) auftreten!). 
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(Die Werte im ersten Oktanten der Abb, 1 erhält man durch Vertauschung der Werte von p, und 9, sowie 


von P(T) durch P+) = 1 — PT) 


Pet Werte des Parameters p(z) bzw. p 
aa Saar 
Rab 


e 
B 
Ei 
u 
ui 
7 
Bi 
7 


BNANAME 


7 n 
Pr TZ Zr 
PU ITAED AD VA 


Ü) 
03 1 _ 
rLrTIr 
e Baal 
Eugen 
= | 
erg wen 
‚es ZanpE® 
[7] 07 02 03 04 05 06 07 08 09 79, 
Bild1l 


Wir vergleichen im folgenden die Hypothesen 
Pp=P, d.h. Po) = 0,5 und 9, > 9, d.h. Pa = 
p > 0,5. (Der Vergleich von P7, = 0,5 und Pn=P 
<0,5 läßt sich analog diskutieren.) p kennzeichnet 
er bestimmte, von p, abhängige Differenzenvertei- 

Zur Beurteilung eines bestimmten Treppenzuges 
konstruiert man nun nach A. WALD in der xy-Ebene 
die Grenzgeraden ax + b mit \ 

In2gq 


ap) = m?» 


(> fürn >- |; A= Ira (3), 
& 
In A ” 
b(p, A) = in2p (> 0 für p> und O0 <a< = 
A (4) 
(s. [4]; Kapitel 6, insbesondere Formeln (6:7) und 


(6:8) für a=ß, wel, MS t=2c-+y, 


| 0,9268 


0,1774 | 0,3267 | 0,4541 | 0,5641 | 0,6600 | 0,7444 
0,1649 | 0,3077 | 0,4324 | 0,5424 | 0,6400 | 0,7273 


vn 
0,1429 | 0,2727 | 0,3913 | 0,5000 | 0,6000 | 0,6923 
0,1379 | 0,2646 | 0,3815 | 0,4896 | 0,5900 | 0,6834 
0,1330 | 0,2566 | 0,3718 | 0,4793 | 0,5800 | 0,6744 
0,1284 | 0,2489 | 0,3623 | 0,4691 0,5700 | 0,6654 
0,1239 | 0,2414 | 0,3529 | 0,4590 | 0,5600 | 0,6563 
0,1196 | 0,2340 | 0,3438 | 0,4490 | 0,5500 | 0,6471 
0,1154 | 0,2269 | 0,3347 | 0,4390 | 0,5400 | 0,6378 
0,1113 | 0,2199 | 0,3258 | 0,4291 0,5300 | 0,6285 
0,4194 | 0,5200 | 0,6190 
0,1037 | 0,2065 | 0,3085 | 0,4096 | 0,5100 | 0,6096 


1 % 1 1 

‚95 0 0,6786 | 0,8261 | 0,8906 9268 

n u) 0,5000 | 0,6920 | 0,7940 | 0,8570 | 0,9000 
0,85 0 0,3864 | 0,5862 | 0,7083 | 0,7907 | 0,8500 
0,8 0 0,3080 | 0,5000 | 0,6320 | 0,7270 0,8000 
0,75 0 0,2500 | 0,4290 | 0,5630 | 0,6670 | 0,7500 
0,7 0 0,2059 | 0,3684 | 0,5000 | 0,6087 | 0,7000 
0,68 0 0,1910 | 0,3469 | 0,4766 | 0,5862 | 0,6800 
0,66 0 

0,64 0 

0,62 0 0,1535 | 0,2897 | 0,4115 | 0,5210 | 0,6200 
0,6 0 
0,59 0) 

0,58 0 

0,57 0 

0,56 0 

0,55 0 

0,54 0 

0,53 0 

0,52 0 0,1074 | 0,2131 | 0,3171 

0,51 0 

0,5 0 0,1 0,2 0,3 0,4 


ER 


mw 


7 


Fe 


| 0,5 0,6 
durch Ersetzung 


a; bzw. r,—= y). Hat dann der Endpunkt des Treppen- 
zugesim Versuchsstadium N speziell die Koordinaten 
x,y, dann ist im Falle ae—b <y<azx-+b der 
Versuch fortzusetzen. 


Für y=ax+b ist wenigstens mit der Wahr- 
scheinlichkeit 1 — « der Wert prichtig (d.h. 9, > 7,). 
Für ax—b > y ist höchstens mit der Wahrscheinlich- 
keit & der Wert p richtig; eine Fortsetzung des Ver- 
suches ist praktisch sinnlos, sofern man hofft, dadurch 
noch p sichern zu können. Weiterhin seien für x und y 
beliebige positive Werte zugelassen. Symmetrische 
Mar (Ps < Pı) ergibt bei festen Werten p und & 

ild 2. 


— P mit mind.1-@gesich.; 
Pz>Pı 


fortsetzen 
PuT-pmit höchst.@gesich.; 
Pr=Pı 


Y 
(7) 


.———— fortsetzen 


*—— T-p mit mind. 1-@ gesich.; 
Pe <Pı 


x?!) 


Bild 2 


‚ Bei festem p (d. h. 9, — p, ist nach Bild 1 als Funk- 
tion von p, gegeben) erhalten wir im Prinzip den Plan 
nach Bild 3. 

Bei festem & (d.h. die Aussagesicherheit (1 —«)- 
100% ist gegeben) dagegen sieht er entsprechend 
Bild 4 aus. 

Wir finden, daß (für y > x) die oberen Grenzgera- 
den sich schneiden. Für Folgen, die in dem schraf- 
fierten Gebiet — etwa in P — enden, haben wir die 
beiden Aussagen: Eine gewisse, durch p=0,8 ge- 
kennzeichnete Differenz ist wenigstens mit 1—« ge- 
sichert, eine (bei gleichem 2,) kleinere, durch 0,7 
gekennzeichnete jedoch nicht. Der zweite Teil der 
Aussage interessiert dann aber i. a. nicht mehr; wir 
können die Strecke g fortlassen. 
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p=0,7 über y=X 
7 ) 2=03 unter y= x 


- 


| @&= 2 entspricht 90% Auss.sich, 


ass entspricht 95% Auss.sich. peg3-- - a3 EIISDRETE 


90% Auss.sich. 
x(2) 
Bild 3 
Diese Feststellung erlaubt die Herstellung eines ein- x(%) 
heitlichen Planes für beliebige « aus 0 <«< 1/2 und Bild 4 
beliebige paus 0<p<1. Dazu ersetzen wir die un- 
übersichtliche Schar der oberen Grenzgeraden durch 
ihre mit 9 bezifferte Einhüllende, von der dann allein Tabelle 2 
die (auch nicht eingezeichneten) rechten Halbtangen- ’ — 
ten interessieren. (Man überzeugt sich leicht davon, l B 1 
daß für die unteren Grenzgeraden entsprechende Ver- ae: s N Pin ler ze ı 
infachungen nicht möglich sind.) 
einfachung g 
Die Einhüllende ergibt sich aus (11) mit (12) (13) 1% 

im Grenzfalle ( =1—.p) | => 
F&yp)=y-nm2p+tzin2g—nA=0 (6) 00 1 0 1,442689 | 1 0 
Se 0,05 | 0,952381| 0,094914| 1,898284 | 0,95 0,052631 

IE go) O1 10,909091 0,233992| 2,339920|0,9 [0,1111 
rt 015 |0,869565| 0,426342| 2.842278 | 0,85 | 0,176471 
Bed}, 0,2 ‚0,833333| 0,687030) 3,435151|0,8 |0,250000 
durch Elimination von (x und y > 0) 0,25 |0,800000| 1,037657 4,150626 | 0,75 | 0,333333 
Eee) 7 0,3 | 0,769231 1,508880 5,029599 | 0,7 | 0,428571 
P=aıy tr 0,3510,740741| 2,145042| 6,128690 | 0,68 |0,470588 
zu 0,4 |0,714286 3,011232 7,528 080 | 0,66 | 0,515151 
vn 42m -mA... (sg) 045 0,689655| 4,207121| 9,349158 | 0,64 |0,562500 
h 2 -+y a) 0,5 | 0,666666 5,386 162 11,772323 | 0,62 | 0,612903 
oder mit 0,55  0,645161 8,299006 | 15,089101 | 0,6 |0,666666 
z=&-.nmA, y=n:mA ....(9) 0,6 |0,625000, 11,873664| 19,789440 | 0,59 | 0,694915 
le io onabhangie vona: 0,65 |0,606061| 17,378053| 26,735466 | 0,58 | 0,724138 
0,7 ,0,588235| 26,300958 37,572798 | 0,57 | 0,754386 
N EN et 0,75 |0,571429| 41,858518 55,811358 | 0,56 | 0,785714 
l+c 0,8 |0,555555| 71,839080 | 89,798851 | 0,55 | 0,818181 
N (10). 0,85 |0,540541| 139,710716| 164,365548 | 0,54 | 0,851851 
en nd mi ] 0,9 |0,526316 | 342,596117 | 380,662352 | 0,53 | 0,886 793 
a er a ae 0,95 |0,512821 | 1464,919038 | 1542020040 | 0,52 | 0.923077 
Die Auswertung ergibt Tabelle 2 und den Plan nach 1,0 10,5 | [6°°) (6) 0,51 | 0,960784 
Tabelle 3 und Bild 5. | 0,5 1 
Tabelle 3 
& = 0,4 (60%) x = 0,3 (70%) « = 0,2 (80%) & = 0,1 (90%) 
c In A = 0,40547 In A = 0,847 In A = 1,38629 In A = 2,19722 
x y & y ® | y ® y 
0,0 ) 0,585 0 1,222 0 | ..2,000 0 |. 8,170 
0,05 0,039 0,770 0,080 1,608 0,132 2,632 0,209 | 4,171 
0,1 0,095 0,949 0,198 1,982 0,324 3,244 0,514 5,141 
0,15 0,173 1,153 0,361 2,407 0,591 3,940 0,937 6,245 
0,2 0,279 1,393 0,582 2,910 0,952 4,762 1,510 7,548 
0,25 0,421 1,683 0,879 3,516 1,439 5,754 2,280 9,120 
0,3 0,612 2,039 1,278 4,260 2,092 6,973 3,315 11,051 
0,35 0,870 2,485 1,817 8,191 2,974 8,496 4,713 13,466 
0,4 1,221 3,052 2,551 6,376 4,174 10,436 6,616 16,541 
0,45 1,706 3,791 3,563 7,919 5,832 12,961 | 9,244 20,542 
0,5 2,387 4,773 4,986 9,971 8,160 16,320 12,933 25,866 
0,55 3,365 6,118 7,029 12,781 11,505 20,918 18,235 33,154 
0,6 4,814 8,024 10,057 16,762 16,460 27,434 26,089 43,482 
0,65 7,046 10,840 14,719 22,645 | 24,091 37,063 38,183 58,744 
0,7 10,664 15,235 22,277 31,824 36,461 52,087 57,789 82,556 
0,75 16,972 22,630 35,454 47,272 58,028 UN AK 2 U bege 91,972 122,630 
0,8 29,129 36,411 60,848 76,060 99,590 124,487 157,846 197,308 
0,85 56,649 66,645 118,335 | 139,218 193,680 | 227,858 | 306,975 | 361,147 
0,9 138,913 |, 154,347 290,179 | 322,421 474,938 527,708 | 752,759 836,399 
0,95 593,981 | 625,243  1240,786 | 1306,091 2030,803 | 2137,687 3218,749 | 3388,157 
1,0 (ee) [6 0) (6°) [e°) © | © [60) 60) 


(Hinter « steht in Klammern jeweils die Aussagesicherheit in %: (1—s) 100%!) 
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j 2 _ Tabelle 3 (Fortsetzung) E 


er & = 0,05 (95%) & = 0,01 (99%) & = 0,0027 (99,739 « = 0,001 (99,99 
R In A = 2,94444 nA= Sons hA= AR nA= a 
3 © | 4Y © Yy © Yy © Y 
en FE a Tr ae EA 
0,0 0 4,248 0 6,629 0 8,529 On 9,957 
0,05 0,280 5,589 0,436 8,723 0,561 11,222 0,656 13,111 
0,1 0,689 6,890 1,075 10,752 1,383 13,833 1,616 16,161 
- 0,15 1,255 8,369 1,959 13,061 2,521 16,803 2,945 | 19,631 
0,2 2,023 | 10,115 3,157 15,785 4,062 - 20,308 4,745 23,726 
0,25 -3,055 | 12,221 4,768 19,073 6,134 24,538 7,167 28,667 
0,3 4,443 | 14,809 6,934 23,112 8,920 29,734 10,422 34,738 
0,35 6,316 | 18,046 9,857 28,162 12,681 36,232 14,815 42,329 
0,4 8,866 | 22,166 | 13,837 34,592 17,802 44,505 | ° 20,798 51,995 _ 
0,45 12,388 | 27,528 | 19,332 42,961 24,872 55,270 29,058 64,572 
0,5 17,332 | 34,663 | 27,048 54,095 34,798 69,596 40,654 81,309 
0,55 24,436 | 44,429 | 38,135 69,336 49,062 89,204 57,319 104,217 
0,6. 34,961 | 58,269 | 54,561 90,935 70,195 116,991 82,008 136,681: 
0,65 51,169 | 78,721 | 79,854 122,853 102,736 158,055 120,026 184,655 
0,7 77,442 110,631 | 120,856 172,652 155,986 | 222,123 | 181,654 259,506 
0,75 123,250 | 164,333 | 192,345 256,460 247,459 | 329,946 | 289,106 385,475 
0,8 211,526 264,407 | 330,109 412,637 424,698 | 530,873 496,175 620,218 
0,85 411,370 | 483,965 | 641,988 755,279 825,942 971,696 964,947 1135,232 
0,9 1008,754 | 1120,838 | 1574,270 | 1749,189 2025,360 |2250,400 2366,226 2629,140 
0,95 | 4313,366 | 4540,385 | 6731,479 | 7085,767 8660,308 | 9116,114 | 10117,829 10650,347 
1,0 [6) 6°) oo 6) [6°) [6°) 0) (6°) 


ER 


(Hinter & steht in Klammern jeweils die Aussagesicherheit in %: (1—o) 100%!) 


Zur Beurteilung eines Treppenzuges entnimmt man 
der durch seinen Endpunkt gehenden Einhüllenden 
die Werte x und p. Die durch p gekennzeichnete Dif- 
ferenz ist dann mindestens mit 1 — « gesichert, sofern 
etwa für den Endpunkt y > x ist (sonst s. Bild 2). 
Man kann auch durch einen hinreichend hoch liegen- 
den Endpunkt die Tangente an eine zu einem ge- 
wählten « gehörende Einhüllende zeichnen und 7 (evtl. 
interpoliert) dem Berührungspunkt entnehmen. 


20 40 P60oX80 100 120 140 160 780 200 
Bild 6 


Mit Rücksicht auf die stochastische Unabhängig- 
keit der Folgen findet man leicht (N=z2+y+P-+M) 


De! x er IE 
ner P+y 


REBSTefe Yy 12 
Meat rer 


CH 
phisch dem Nomogramm (Bild 6) entnehmen. Es ent- 
hält p als Parameter an den Radialstrahlen, und in 


Die Ausdrücke der Form p = = T — lassen sich gra- 


ihm lassen sich auch in der angedeuteten Weise die 
Mittelungsprozesse ausführen. So ermittle man nähe- 
rungsweise p,. Zusammen mit pkann man dann dem 
Bild 1 die wenigstens mit der Aussagesicherheit 
1 — « absicherbare Differenz p,—p, entnehmen. 
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On the Displacement-Function for the Axially 
Symmetrical Visco-Elastic Vibration 


The displacement-function for the solid of revo- 
lution subjected to the axially symmetrical vibration 
has been found by me [1] and M. PREDELEANU [2]. 
The present article gives its extension to the case of 
the visco-elastic vibration, provided the stress-strain 
relation is governed by Voıgr’s law. 

This law is written 
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and the displacement-equations become 
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u and ıw being displacement-components in direetions 
to r and 2 axes. 

Now it will be eonvenient to introduce the two 
notations 


0 0° aa 
a Er oe dt 


8)» 


EN 


a and b being termed wave-slownesses by HAMILTON. 


Assuming a displacement-function x for. the system 
under consideration, we may put 


,- 1 o\ 0% 
u, [14r5)0% 

1 & ‚0\o 
=, a-n {14% 327 


where x and x’ are constants to be determined. -Sub- 
stituting equations (5) into the first of equations (2) 
and rearranging the result, we have 


a Mh 
2uorxöt\!-+u 


8); 


EDEL ZEN, 


’ ! x , ö 0° 
+0) [a4 „+ (A 4 2 
X"H+w\ 0° 
7 «3 e - anal, 
This equation is identically satisfied by taking 
Es ; 
== ‚Bsp: re: 


In virtue of equation (6), the second of equations (2) 
can be rearranged, with equations (5), to the form 


Here ern: 


the equation for the displacement-function Hr 
The displacements then become 


1 X“+w 0\ 8% 
=—_—|l 
% Er FR =) Or 02’ i 
(8), 
en [ XH+u 8\ 8 
u ,, 1 „Di nr 


and the corresponding stresses become 


Ber Yo\, « w2 
Mu AN a 
ö X 6 .T Er) 
BE n 
x (14 5) al u 
IR ah. X & $ wWö 
2(1—») {1 
Y+u ) 


x 


u 


ll A 
tu lol * 


The dilatation and rotation take the forms 


1—2v ö 


R 1—v & 
A 


01%» 
so that, with equation (7), 


D;A4A=0, 


If, in equations (7), (8), and (9), we put#’ = u’ =, 
then these equations reduce to those derived by me [1] 
and PREDELEANT [2]. 

It is added that the calculation begins with the 
equations 


KO=R. 


or? 2, Ir U 

— ir 9e 
a ' & * 7 , er 
oz 2: m _ MW, 


U, W being some functions of y, z and t, and the calcu- 
lation procedure is the same as that given in Love’s 
Elasticity [3]. 

In the case of thermo-elasticity governed by 
DuHamer’s law, similar calculation will afford 


Vy=— pe Ira, 
en 177935 
24 Or & 


w =;- un Tzi— zaler2 te | Tiazz 


T being a temperature distribution, E_ Young’s 
modulus, and c thelinear coefficient of thermal expan- 
sion of the solid. 
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Über die Abrundungsfehler bei der Iteration für 
y=ye. 
Das Newroxsche Iterationsverfahren zur Bestim- 


mung von Nullstellen liefert bekanntlich, angewandt 
auf die Funktion f(y) = y* — x, für 
ya  @>0;n=23,34..).: 


eine Folge von Näherungswerten y; gemäß der Vor- 
schrift 


Im: 
Yk+1 lg (n —1) %Yk (k = 0, 1, 2,» .) 


Yr 


Wählt man y, > y, so gilt: 
Ra DEF IR ay ..(B). 


Dabeisind &,...,%%...,%1i. a. nicht abbrechende 
Dezimalzahlen; mit einer Rechenmaschine können 
aber, bedingt durch die Kapazität des Rechenwerkes 
[2 p Stellen] und der Speicher [p Stellen], lediglich 
solche Zahlen verarbeitet werden, die von der (o p+1)- 
ten Stelle an (o=1,2,...,r) Nullen aufweisen.!) 
Die theoretische Iterationsvorschrift (2) ist also 
allein nicht ausreichend, sie muß vielmehr durch 
Abrundungsvorschriften — und zwar vor allem für 
die Divisionen?) (x — q‘®) 


Pina, My”, 
es ER Sy gr=n 


(4) 


— so ergänzt werden, daß man als Ergebnis der 
Rechnung eine nicht zu verbessernde p-ziffrige 
Näherung für y erhält. 

Dazu ist zunächst zu präzisieren: 


x sei eine 9-ziffrige, positive Zahl in der Darstellung 
(0, pP): 
2 m 10727 min — 1,2,..., 107 7%) 


gesucht ist, ebenfalls in der Darstellung (0, p), die- 
jenige p-ziffrige, positive Zahl y, für die 


meyjz zus tr... .(6) 
ausfällt. 

Ferner wird — in Verallgemeinerung von (6) — zu 
jeder positiven Zahl 2 (mit O0<2z<1) durch Ab- 
rundung nach unten eine o p-ziffrige Zahl [2],» de- 
finiert: 

[z2]op entsteht aus 2, indem von der (op + 1)-ten 
Stelle an jede Ziffer durch Null ersetzt wird; es gilt 
also 

Bra 2er 10 Pe, (7). 
Jetzt läßt sich zeigen: Geht man von einer p-ziffrigen, 
positiven Zahl %, > y aus°®), so ergibt sich y” als 
letztes Element yx einer endlichen, monoton fallenden 
Folge »-ziffriger, positiver Zahlen %%, in der Dar- 
stellung (0, 2): 


H>Yyı> ">y>Yrı>'>Yy-ı >yZYrd). 


Man erhält diese Folge, wenn das Newtonsche Itera- 
tionsverfahren (2) in folgender Weise‘) modifiziert 
wird (2 = 99): 


go == [ge-» n le (o _— 17 2, auele 


Yrrı = [9 + (n — 1) Yr) : Rp 


Die Iteration wird fortgesetzt, solange % > Yk+ı ist; sie ist beendet, und zwar 
mit dem Wert y, = Yk = y,, sobald diese Bedingung nicht mehr erfüllt ist. 


ı) Man wird unter allen Umständen bemüht sein, r und damit den 
Rechenaufwand so klein wie möglich zu halten. 

2) Der Übersichtlichkeit halber wird bei den Quotienten g(0) ein 
besonderer Index %k weggelassen. 

3) Auch % = 1 ist möglich und zumeist sinnvoll. 

4) Die Iterationsvorschrift (9) gilt unverändert auch für Dual- 
zahlen. 


a) 


Beweis. Es sei yp+1 der theoretisch, d.h. gemäß (2) 
aus % > y sich ergebende Wert, dann gilt: 


N 
ne 43% 
Rn >0E TE RR (1 


N ee] ——. -1fzn-o n-o 
Yk+ı — Y= En 2% {7% —Yy }>0 (10.2) 
> GE 
also 
Yr > Ycsı >Y (k=.0,12;, 2.97, Ra 
Aus den Abrundungsvorschriften (9.1) schließt man 


ge» 10 Rr-e+1)» 


see ia 
Yk Yk 
(e =1,2,2..,n—1) 


(12.1) 


und kann nach Division dieser Ungleichungen durch 
% °!(o=1,2,...,n—1) über og summieren: 


(0) ; 
< aa Zgn-2 E10, rn (12.2). 
k 3 


Da (9.2) besagt, daß 
08 (ud +1) 5) in — er SU 10 
(12.3) 


ist, erhält man schließlich durch Zusammenfassung 
von (12.2) und (12.3) 


0 


IN 


1 2 = ® 2 
7 gatmls ea ee 
(12.4). 


Es ist also %%+1ı = [Yr+ı]p; bei jedem Iterations- 
schritt wird insgesamt um weniger als eine Einheit 
der letzten Stelle nach unten gerundet. Wegen (11) 
ist andererseits Y,,ı um mindestens 10-? kleiner als 
Yr; die Folge %7, ist also endlich, es gilt (8) und ins- 
besondere yx_ı >yk >y 2 Y%&d.h.yk = y. (Aus 
Yk < y folgt entsprechend Y% <y<yx,;ı und damit 
YK < YK+ı, also nicht mehr % > Yk+ı-) 

Beispiele zeigen, daß sich bei der üblichen Ab- 
rundung, d.h. wenn stets in der p-ten Stelle anstatt 
gemäß (9) gerundet wird, häufig schlechtere Näherun- 
gen ergeben. Für n = 3 erhält man (%, = 1;p = 10): 


es Rundung in der 


© Y zehnten Stelle 


0,75000 00000 
0,00800 00000 
0,00001 00000 
0,00000 00001 


0,90856 02964 
0,20000 00000 
0,02154 43469 
0,00046 41588 


0,90856 02963 
0,19999 99998 
0,02154 43455 
0,00046 41120 


. (9.1) 
. (9.2) 


. (9.3) 


Verfasser: Dr. JULIUS ALBRECHT, Institut für An- 
gewandte Mathematik der Universität 
Hamburg, Rothenbaumchaussee 67/69 
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Fehlerabschätzungen für das Iterationsver- 
fahren von Schulz zur Bestimmung der In- 
versen einer Matrix 


Ein Iterationsverfahren zur Bestimmung der In- 
versen einer Matrix ist von G. ScHuLz!) angegeben 
worden, für das R. AnsoRGE eine Fehlerabschätzung 
bewiesen hat?). Die Abschätzung von ANSORGE hat 
den Nachteil, daß sich der zur Durchführung eines 
Iterationsschrittes benötigte Rechenaufwand etwa 
um die Hälfte vergrößert, sobald man den Ite- 
rationsfehler abschätzen will, wenn die zu inver- 
tierende Matrix und die Matrizen der Näherungsfolge 
nicht symmetrisch sind. Man kann jedoch völlig 
zwangsläufig zu Fehlerabschätzungen kommen, denen 
dieser Nachteil nicht anhaftet, wenn man zum Beweis 
die funktionalanalytische Beschreibungsweise ver- 
wendet. Da dem Iterationsverfahren von SCHULZ trotz 
seiner Wirksamkeit bis jetzt noch nicht die ihm ge- 
bührende Beachtung in der numerischen Mathematik 
geschenkt worden ist, sei es gestattet, an dieser Stelle 
nochmals auf das Verfahren einzugehen. 


1. Bemerkungen zur Durchführung 

des Iterationsverfahrens 

Die Inverse einer nicht singulären, n-reihigen, 
quadratischen Matrix C ist Lösung der Matrizen- 
gleichung 


OR SHE ee (1) 
und damit auch der Gleichung 
ZEIT OD SOHN re (2). 


Beim Iterationsverfahren von SCHULZ wird aus- 
gehend von einer nicht singulären Matrix X, eine 
Folge von Näherungsmatrizen nach der Vorschrift 


Xp = Xn (2 E—COX,) - (8) 


gebildet. Sämtliche Matrizen dieser Folge setzen wir 
als nicht singulär voraus. Die Rechenarbeit bei jedem 
Iterationsschritt besteht im wesentlichen aus zwei 
Matrizenmultiplikationen. Verwendet man zur Ma- 
trizenmultiplikation die von 8. FALk?) vorgeschlagene 
. Anordnung, so kann der Iterationsprozeß nach dem 
folgenden Rechenschema ausgeführt werden. Das 
Schema ist für zwei Iterationsschritte angegeben 
ng wobei der erste Schritt stark umrandet wor: 
en ist. 


Die Kontrolle der Rechnung kann durch Zeilen- 
oder Spaltensummenkontrollen erfolgen. 
5} ’@ = . 
2. Konvergenz des Iterationsverfahrens 


Zum Beweis der Konvergenz der Iterationsfolge 
gegen die Inverse der Matrix © berechnen wir 
E—COXn+1. Aus (3) ergibt sich: 

E— OXn +1 =E— OXn (2 E— CXn) 
= E—20X, +0X,0X, 
= (BC 22222. li), 


!) G. SCHULZ, Iterative Berechnung der rezi { ix, 2 
2% (ss), A g eziproken Matrix, ZAMM 
) R. ANSORGE, Über ein Iterationsverfahren von G, S 
| ® . SCHULZ 
aus der Reziproken einer Matrix. ZAMM 39 (1959), S. 164 bis 


2 

>) S. FALK, Ein übersichtliches Schema für die Matri - 
plikation. ZAMM 31 (1951), 8. 152—153, Siehe auch Ne 
Matrizen, 2. Aufl., Berlin, Göttingen/Heidelberg 1958, Springer-Verlag 


S. 17. 
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Diese Gleichung zeigt die quadratische Konvergenz 
des Verfahrens. Aus (4) kann sofort 3 


(E— CXn+ı) = (BC) 


gefolgert werden. 

Wir führen nun’im Raum der n-reihigen, quadrati- - 
schen Matrizen eine Norm ein. Unter der Norm einer 
Matrix A = (a;j) wollen wir hier die maximale Zeilen- 
summe der Beträge der Elemente verstehen: 


1 


ae) 


n 
All = max 3 ur. 
ik — 


Auch die Verwendung eines anderen Normbegriffes 
kann gelegentlich von Nutzen sein. 


Setzen wir 
g=1IE—OX ll -» - +. (6) 


so erhalten wir unter der Voraussetzung 
ee a 


NE— CXrsıll>0 fürn>o ... (8). 
Da für zwei Matrizen des Raumes 

lAa2lls|Al-|Bl - - - 
ist, können wir 


Io — Xa+ı]| = ||0* (E— OXn+1)]|| 
s Je ]|- ||E— CXn+a]| 


folgern, womit wegen (8) unter der Voraussetzung 
q<1 die Konvergenz der Iterationsfolge gegen die 
gesuchte Inverse bewiesen istt). 


Damit die Voraussetzung q < 1 erfüllt ist, muß die 
Ausgangsmatrix X, bereits eine „gewisse Näherung“‘ 
für die Inverse ©! sein. Das dürfte auch der Grund 
sein, weshalb diesem Verfahren bis jetzt nicht die 
genügende Aufmerksamkeit geschenkt wurde. In 
vielen Fällen, z.B. wenn die Matrix C überwiegende 
Hauptdiagonalelemente besitzt, wird es nicht schwie- 
rig sein, eine Matrix X, anzugeben, mit der die Voraus- 
setzung (7) erfüllt ist. Man kann sich auch durch 
grobes Rechnen nach einem anderen Verfahren eine 
Ausgangsmatrix X, verschaffen und diese Matrix 
dann in der angegebenen Weise iterativ verbessern. 
Der Vorteil des Iterationsverfahrens von ScHuLz liegt 
in der guten Konvergenz und der Möglichkeit einer 
Fehlerabschätzung. 


aus (5) 


510) 


3. Fehlerabschätzungen 
Aus den beiden Gleichungen (2), (3) finden wir sofort 
X— Xn4ı =2 X — 2X, — XCX + X,C0Xn 
= (X — X.) + (X — AXCX) — X, + Xn COX. 


Berücksichtigen wir, daß X als Lösung der Glei- 
chung (1) auch der Gleichung 


20=E 
genügt und damit auch 
XCX, = A 


gilt, so ergibt sich weiter 
X — Inıı =(X— X.) — XCX, + X, 0X, 
= (X — Xn+ı) + (An+ı — Xn) 
— (X — Xn+1) COX 
u ORR En An) COX, 
(X An+ı) OXn == (Xn+ı — x) (E me OX,) « 


‘) Eine notwendige und hinreichende Bedingung für dieK . 
genz der Iterationsfolge findet man bei G, es (8. Fußnote 1). 


Diese Bedingung i ) i 
Bodantrz gung ist jedoch vom numerischen Standpunkt kaum von 


Ele Suä Ba BE N La a a ucıa Zul ann u Eu m um ud na nn aan 
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Unter den gemachten Voraussetzungen existiert die 
Inverse von COX, so daß wir schreiben können 


X — Xn+ı = (Xn+ı — X) (E— OX,) (COX)! 
Durch Übergang zu den Normbeträgen erhalten wir 


mit (9): 
R X — An+ıl| 
-S |ICXm= ||- |B— CXn||- | nt —&u|| _ - (10). 


Nun läßt sich bekanntlich (E — A)-! analog zur geo- 
metrischen Reihe gewöhnlicher Zahlen in eine Matri- 
zenreihe, die sogenannte NEuUMANNsche Reihe, ent- 
wickeln 


EA ER I Arne 


die für ||A|| < 1 konvergiert. Daher gilt mit ||2]] = 1 
unter Beachtung von (9) die Abschätzung 


1 
1— Ill 
Mit Hilfe von (11) wollen wir die Norm von (OX,) 1 
abschätzen. Es ist 
(OXn)ı={E—(E— COX}, 
und unter der Voraussetzung (7) folgt aus (5) 
|#— CXn|<1, 


a — 3]|< Sal Sllalt = an. 


so daß (11) anwendbar ist. Wir erhalten somit für 


IE — OXn|| r a 5 
Is ee Xn||| (12). 


Dafür können wir wegen (5), (6) auch schreiben: 


ie ie 


Will man sich die Berechnung von An+1 — An 
ersparen, so kann man ihre Norm noch abschätzen. 


(13). 


"Aus (3) finden wir 


Anti — An = An (E— CR, 
In — || S ||| | — CA < "||. 


Für die Abschätzungsformeln (12), (13) ergibt sich 
dann: 


E—CX,|? ® 
IX—-Hlls; ran |||] ao 
® fen 
IX —&r+ıll s; Anl ..(18). 


Für das Iterationsverfahren von ScHuLz können 
diese Fehlerabschätzungen ohne große Mehrarbeit 


‘durchgeführt werden, da die Berechnung von g bzw. 


|? — OX„|| kaum zusätzliche Rechenarbeit erfordert; 
denn die Matrizen HE — COX, werden faktisch durch 
das Verfahren geliefert, und es-muß nur noch die 


Norm dieser Matrizen sowie die Norm der Iterations- 


matrizen bzw. gewisser Differenzen von Iterations- 
matrizen bestimmt werden. 


4. Beispiel und Vergleich der Fehler- 
abschätzungen. 


Zur Demonstration soll das Iterationsverfahren von 
SCHULZ für eine dreireihige, quadratische Matrix 
durchgeführt werden. Es werden zwei Iterations- 
schritte berechnet. Das Rechenschema ist nach den 
im ersten Abschnitt dieser Arbeit gemachten Be- 
merkungen von selbst verständlich. Die kursiv ge- 
druckten Spalten enthalten die Zeilensumme der Be- 
träge der Elemente der Matrizen X, bzw. E— OX,; 
die maximale Zeilensumme, welche die Norm der 
Matrix bestimmt, ist zusätzlich fett gedruckt. 


Dem Rechenschema entnehmen wir folgende für die 
Fehlerabschätzungen benötigte Zahlenwerte: 
Ill = 6,7 IGl| = 691 
q= ||E— CAX,|=09<1ı |E— CXj| = 0,51. 
Leicht lassen sich die Gleichungen 
||Xı — X,|| = 0,49 IX — X,|| = 0,0977 


bestätigen. 

Die Ergebnisse der Fehlerabschätzungen stellen 
wir in einer Tabelle zusammen und vergleichen sie 
mit dem exakten Fehler: 


IX — ll | 1IX— Zoll 
Abschätzung nach (12) 4.41 0,102 
Abschätzung nach (13) z 0,417 
Abschätzung nach (14) 54.27 3,670 
Abschätzung nach (15) ai 23,87 
Exakter Fehler 0,15 0,0525 


Für den ersten Iterationsschritt fallen naturgemäß 
die Abschätzungen (12), (13) sowie (14), (15) zusam- 
men. Die Abschätzungen nach dem ersten Iterations- 
schritt sind sehr schlecht. Das liegt daran, daß von 
einer Matrix X, ausgegangen wurde, die zu einem re- 
lativ großen Zahlwert von q führt. Das gewählte 
X, ist eine numerisch nicht günstige Ausgangsmatrix. 
Das Beispiel wurde aber so angelegt, um die Unter- 
schiede zwischen den einzelnen Fehlerabschätzungen 
deutlicher hervortreten zu lassen. 


Beim zweiten Iterationsschritt liefert die Fehler- 
formel (12) bereits sehr gute Ergebnisse; der Zahl- 
wert für |E— CX}|| ist mit 0,51 auch numerisch 
günstig. In Formel (13) geht bereits der numerisch 
ungünstige Wert von g ein, und für die Abschätzun- 
gen (14) und (15) verschlechtern sich die Ergebnisse 
aus diesem Grunde noch mehr. Hätten wir X, zur 
Ausgangsmatrix gemacht, so hätten wir für die Ab- 
schätzungen (12), (13) den Wert 0,102 und für (14), 
(15) den Wert 3,670 erhalten. 


Die Abschätzung (15) ist zweifellos numerisch am 
einfachsten, dafür aber auch am schlechtesten; die 
Anwendung von (14) erfordert jedoch kaum mehr 
Rechenarbeit, man hat lediglich bei den weiteren 
Iterationsschritten die Norm einer durch das Verfah- 
ren faktisch gelieferten Matrix zusätzlich zu bestim- 
men. Entsprechendes gilt für die Formeln (12) und 
(13). Bei der Abschätzung (12) hat man zusätzlich 


Bon seelemä: | 11 0200 | 0,8 09 —,10 004 |0,1| 
E00 2.0.2%1.0,9 0065 136 0,10 |051| 
| 9 01 21 |02 [000 —0,05 1,08 10,09 | 
| 27 1-09 00 18|--0,99 0,00 1,98 | 2:97 | —0.9999 0,0000 1,9998 
67 | 37 10-—2.0| 397 094 —2,0016,91 | 3,9973 0,9814 —1,9952 
150 | 28 21-11] 297 0,89 —0,99 | #35 | 2.9947 0.9629 —0,989 
Tea 09702. 001 1.01, =0.10770,08 | 
a NT 16. 02.1. 008: —0,64 0,10 
ı 1-1| 00 0,1 —0,9 | 0,01 —0,05 —0,97 | 


gegenüber (14) noch 


die Matrix Xn+1 — z zu 


rechnen, was aber durch die günstigeren Resuli tate 


der Abschätzung i. a. zu empfehlen ist; das gilt auch 


für die Formeln (13) und (15). Man kann sagen, daß 
die Abschätzung nach (12) immer vorteilhaft sein 
wird, wenn q nicht sehr klein ist; in diesem Falle 
unterscheiden sich die angegebenen Fehlerabschätzun- 
gen i. a. nur sehr wenig voneinander. 


5. Schlußbetrachtungen 

Zum Schluß sollen die angegebenen Fehlerab- 
schätzungen mit der Abschätzung von ANSORGE ver- 
glichen werden. In der hier benutzten Schreibweise 


lautet die Abschätzung von ANSORGE mit X = (zij), 
Xu= (ef) 

(HR +1) <_Ir+l _ (n+1) _ „in 
Wr at, \=7 Sn a =7|» 


BUCHBESPRECHUNGEN 


Grimsehl, Lehrbuch der Physik.Zweiter Band: 
Elektromagnetisches Feld. 14. Aufl. Herausgegeben 
von Prof. Dr.W.Schallreuter (DirektordesPhysikal. 
Instituts der Universität Greifswald). VIII + 563 S. 
m. 680 Abb. Leipzig 1959. B.G. Teubner Verlags- 
gesellschaft. Preis geb. 16,20 DM. 


Mit der 14. Auflage hat der II. Band des Grimsehl 
eine völlige Umarbeitung erfahren. Man bemerkt mit 
Freude, daß sein Charakter erhalten blieb, seine An- 


'schaulichkeit, das liebevolle Eingehen auf experimen- 


telle Einzelheiten und die vielen Anwendungen der 
Physik in der Technik. Umstellungen und Kürzungen 
ermöglichten es, den Umfang des Bandes trotz rei- 
cheren Inhalts beizubehalten. Neu sind z. B. die Ab- 
schnitte über die molekulare Struktur der Elektrolyte, 
Frequenzwandlung und Filterung, Modulation, Er- 
zeugung und Übertragung ultrakurzer Wellen. Der 
Abschnitt über Elektronenröhren ist erweitert, Tran- 
sistoren werden kurz erläutert. An anderen Stellen 
wurden neue Beispiele eingefügt. Die mathematischen 
Anforderungen an den Leser sind die gleichen geblie- 
ben. Insgesamt hat die Darstellung an Klarheit ge- 
wonnen. Man könnte dieses Lehrbuch den Studenten 
der Technischen Hochschulen uneingeschränkt emp- 
fehlen, wenn das MKSA-System benutzt würde. 
Der Herausgeber hat sich jedoch zur Verwendung des 
Gaußschen Maßsystems entschlossen. Dies empfindet 
der Referent als einen beträchtlichen Nachteil. 


Nachdem 1948 das MKSA-System international 
angenommen war und schließlich 1958 in der DDR 
für verbindlich erklärt wurde, konnte man hoffen, daß 
dem Studierenden der Anfangssemester die Mühe des 
Unterscheidens mehrerer Maßsysteme in Zukunft er- 
spart bleibt. Der vorliegende Band enttäuscht diese 
Hoffnung ein wenig. Es ist zu bedauern, daß diese 
Schatzkammer experimenteller Erfahrung nicht diesen 
bequemeren Zugang erhalten hat. Eine weit größere 
Zahl angehender Ingenieure würde zu diesem sonst 
vorzüglichen Buch greifen. Schließlich sei noch be- 
merkt, daß das hier verwendete Kurzzeichen für die 
Einheit der elektrischen Ladung sowie die Bezeich- 
nung mp. für 10 m nicht der gesetzlichen Festlegung 
entsprechen. 


Dresden H. ZıimmER 


Loo-Keng Hua, Die Abschätzung von Expo- 
nentialsummen und ihre Anwendung in der 
Zahlentheorie (Enzyklopädie der mathematischen 
Wissenschaften, Band ];, Heft 13, Teil I. Heraus- 
gegeben im Auftrage der Akademien der Wissen- 
schaften zu Berlin, Göttingen, Heidelberg, Leipzig, 
München und Wien sowie unter Mitwirkung zahl- 
reicher Fachgenossen.) 1238. Leipzig 1958. B. G. 
Teubner Verlagsgesellschaft. Preis brosch. 13,— DM. 
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bis (15) nicht an. 


Verfasser: Dr. W. Dück, Berlin-Johannisthal, ne Fe 
Ecksteinweg 8 


Die Abschätzung von Exponentialsummen ist eine 
analytische Methode der additiven Zahlentheorie. Als 
solche wird sie im Anschluß an eine Übersicht über 
die elementaren Methoden dargestellt. Anschließend 


‚werden an zahlentheoretischen Problemen haupt- 


sächlich die Primzahlverteilung, das Waringsche, das 

Goldbachsche Problem und die Gleichverteilung be- 
handelt. Der Schwerpunkt liegt auf der modernen 

Entwicklung. 

Dem Verfasser ist es wohlgelungen, die Vielfalt des 
anfallenden Stoffes so darzustellen, daß der prinzi- 
pielle Gehalt und die große Linie in der Fülle des Ein- 
zelnen nicht untergeht, sondern allenthalben klar 
herausgearbeitet und mit dem historischen Verlauf 
verbunden wird. Durch die am Schluß hinzugefügte 
tabellarische Übersicht der wichtigsten in dem Artikel 
behandelten Probleme mit ihren neuesten Ergeb- 
nissen wird dem Leser eine bequeme Orientierung er- 
möglicht. 

Eine Frage bliebe zu stellen: Ist es nicht eine will- 
kürliche Verengung des Begriffs der additiven Zahlen- 
theorie, wenn einleitend bemerkt wird, die Geschichte 
der additiven Zahlentheorie beginne im 18. Jahrhun- 
dert mit den Problemen von Goldbach und Waring? 

Potsdam M. DRAEGER 


Dr. H. Kaufmann, Dynamische Vorgänge in 
linearen Systemen der Nachrichten- und Re- 
gelungstechnik. 2118. m. 105 Abb. und einem 
Tabellenanhang. München 1959. R. Oldenbourg 
Verlag. Preis geb. 26,50 DM. 


Bei der Behandlung linearer Probleme der Rege- 
lungs- und Steuerungstechnik sucht man im Zeitbe- 
reich den zeitlichen Vorgang der Übertragung, im Fre- 
quenzbereich die zugehörige spektrale Verteilung. Mit 
Hilfe der Gewichtsfunktion stellt der Verf. in den er- 
sten 4 Kapiteln die Zusammenhänge zwischen den bei- 
den Beschreibungsarten ausführlich dar und skizziert 
auch deren Anwendungsbereiche. Das Kapitel 5 bringt 
die am häufigsten benutzten Stabilitätskriterien von 
Hurwitz, Michailow-Leonhard und Nyquist. 
Im folgenden Kapitel werden Approximationsver- 
fahren zur Ermittlung des dynamischen Verhaltens 
von Netzwerken angegeben, wobei neben häufig be- 
handelten Methoden die Berechnung von Übergangs- 
funktionen durch Entwicklung der Erregung nach 
Laguerre-Funktionen erörtert wird. Nach Darlegun- 
gen über Kennwerte von Übergangsfunktionen wer- 
den im letzten Kapitel Systeme mit Abtastung be- 
trachtet, die im Zusammenhang mit digitalen Syste- 
men erheblich an Bedeutung gewonnen haben. 

‚Das Buch ist klargeschrieben und gibt einen schönen 
Einblick in die behandelten Probleme, so daß es allen 
denen warm zu empfehlen ist, die auf dem Gebiet der 
Nachrichten- und Regelungstechnik tätig sind. 


Dresden H. KınpLEer 


Buchbesprechungen 


 . Dr. phil. Ludwig Föppl (em. o. Prof. für Techn. 
' Mechanik, Techn. Hochsch. München. Mitgl. d. Akad. 
‘d. Wiss. zu Madrid und München), Elementare Me- 
chanik vom höheren Standpunkt. 174 S m. 
71 Abb. u. 8 Bildtafeln. München 1959. R. Olden- 
bourg. Preis geb. 20,— DM. 


Es ist immer reizvoll, von hoher Warte aus eine 
Rückschau auf wichtige Stationen der Entwicklung 
einer Wissenschaft und ihrer Begriffswelt zu veran- 
stalten und dabei auch die Persönlichkeiten lebendig 
werden zu lassen, die diese Entwicklung maßgeblich 
beeinflußt haben. Der Verfasser, der sich diese Auf- 
gabe gestellt hat, hat sie in liebenswürdiger Weise ge- 
löst. Grundlegende Probleme der elementaren Me- 
chanik, angefangen bei den am Hebel und an der 
schiefen Ebene wirksamen Gesetzen über die Pendel- 
schwingungen, das Gravitationsgesetz und die Pla- 
netenbewegungen, die Vorgänge beim Stoß und am 
Kreiselundandere mehr dienen ihm dazu, aufzuzeigen, 
auf welchem Stand sich die Mechanik zur Zeit eines 
Archimedes,Stevin,Galilei,Kepler,Huyghens, 
Newton, Leibniz, Euler, Bernoulli, d’Alembert 
und Gauß befunden hat und welchen entscheidenden 
Einfluß diese großen Männer bei der Klärung und bei 
der Bildung der mechanischen Grundbegriffe ausgeübt 
und was sie dazu beigetragen haben. Das Buch führt, 
den Stil einer historischen Betrachtung allmählich 
verlassend, bis zu den Differential- und Integral-, den 
Minimal- und Variationsprinzipien der Mechanik und 
gibt zum Schluß auch einen Einblick in die Proble- 
matik, die zur Einsteinschen Relativitätstheorie 
geführt hat, und erläutert kurz ihre Grundgedanken. 


Nicht befreunden kann sich der Referent damit, 
auf welche Art und Weise auf S. 137 und an weiteren 
entsprechenden Stellen der Extremwert einer Funk- 
tion (bzw. des bei der Variationsrechnung auftreten- 
den Funktionals) definiert wird, wonach ‚‚wenn es sich 
um einen maximalen oder minimalen Wert... han- 
delt, für zwei benachbarte Werte der Abszisse x die 
Funktionswerte dieselben‘ sind. 


Dieser Einwand soll jedoch den Wert dieses schönen 
Buches nicht herabsetzen, das vor dem Leser ein in- 
teressantes, durch die vielen ausführlichen Angaben 
aus dem Leben der großen Forscher auch in mensch- 
licher Hinsicht sehr lebendig gestaltetes Kapitel aus 
unserer Geistesgeschichte abrollen läßt. 


Dresden H. HeEInRIıcH 


W. Hahn, Theorieund Anwendung der direk- 
ten Methode von Ljapunov (Ergebnisse der 
Mathematik und ihrer Grenzgebiete, Neue Folge, 
Heft 22). VIII + 142 S. Berlin/Göttingen/Heidelberg 
1959. Springer-Verlag. Preis geh. 28,— DM. 


Die Ursprünge der Stabilitätstheorie liegen weit 
zurück, als Begründer darf man wohl Poincare und 
vor allem Ljapunow ansehen. Obwohl beide schon 
vor mehr als einem halben Jahrhundert grundlegende 
Stabilitätssätze aufgestellt und Untersuchungsme- 
thoden entwickelt haben, fanden sie doch nur einen 
geringen Widerhall. Erst seit zwei Jahrzehnten wer- 
den ihre Gedankengänge (hauptsächlich von sowjeti- 
schen Mathematikern) fortgesetzt. Allerdings wird 
diese Fortsetzung mit ungeheurem Eifer betrieben; 
sie hat zu einer gewaltigen Ausweitung der Theorie 
und natürlich zu einer Vielzahl von Veröffentlichun- 
gen geführt. Die Impulse für das große Interesse, das 
neuerdings Stabilitätsfragen entgegengebracht wird, 
gehen von der Praxis aus: Die Automatisierung, die 
die moderne Technik kennzeichnet, verlangt die Kon- 
struktion vielfältiger Typen von schnell wirkenden 
Steuerungs- und Regelungssystemen. Daß bei diesen 
die Stabilität weit mehr als sonst eine Rolle spielt, ist 
einleuchtend. Nachdem nun die Ljapunowsche 
Stabilitätstheorie ausgebaut und in gewissem Sinne 
abgeschlossen worden ist, war es dringend erforder- 


lich, die vielen Einzelresultate zu sichten und unter 
einheitlıchem Gesichtspunkt zu beurteilen und zu- 
sammenzufassen. Dieser wichtigen Aufgabe hat sich 
der Verfasser (der wie kein anderer dazu berufen war) 
in dem vorliegenden Ergebnis-Band angenommen. 
Sein Verdienst besteht auch darin, große Teile der 
(meist russisch geschriebenen) Originalliteratur erst- 
malig einem breiten deutschen Leserkreis zugänglich 
gemacht zu haben. Der Bericht befaßt sich beinahe 
ausschließlich mit der sogenannten direkten Methode 
von Ljapunow, jedoch nur im Zusammenhang mit 
Stabilitätsaufgaben. Neuere Untersuchungen (z.B. 
von Yoshizawa) haben gezeigt, daß sich diese 
Methode auch noch zu anderen Zwecken verwenden 
läßt; dementsprechend kann man das Stabilitäts- 
problem in einen umfassenderen Problemkreis (der 
ganz allgemein das asymptotische Verhalten der Lö- 
sungen von Differentialgleichungen betrifft) einord- 
nen. Wer sich aber auf das Problem der Stabilität be- 
schränken will, findet in dem Bericht des Verfassers 
eine nahezu vollständige und trotzdem nicht zu weit- 
läufige Darstellung der modernen Forschung und 
ihrer Resultate. In Kapitel I werden die Grundbe- 
griffe behandelt; hierzu gehören die üblichen Stabili- 
tätsdefinitionen. Kapitel IT bringt kurz die Haupt- 
sätze der Stabilität und Instabilität und Kapitel III 
ihre Anwendungs -auf konkrete Probleme. Hierbei 
handelt es sich nicht um Aufgaben aus der Praxis; 
denn leider liegen in dieser Richtung bisher kaum 
Resultate vor. Kapitel IV widmet sich dem wichtigen 
Problem der Umkehrung der Hauptsätze, das darauf 
hinzielt, die Existenz Ljapunowscher Funktionen 
als notwendige und hinreichende Stabilitätsbedingung 
darzustellen. In Kapitel V wird das Wachstum der 
Lösungen in besonderen Fällen genauer bestimmt; 
der Verfasser erläutert u. a. den Begriff der exponen- 
tiellen Stabilität und die Bedeutung der charakteristi- 
schen Zahl. In Kapitel VI wird die Stabilität nach der 
ersten Näherung und die totale Stabilität (Stabilität 
bei dauernd wirkenden Störungen) erörtert. Kapitel 
VII wendet sich dann den kritischen Fällen zu, in 
denen die erste Näherung zur Stabilitätsanalyse nicht 
ausreicht. Kapitel VIII bringt schließlich den Begriff 
der Stabilität in einem endlichen Intervall und andere 
Verallgemeinerungen (Übergang vom Euklidischen 
Phasenraum zu allgemeineren Räumen, Stabilität bei 
partiellen Differentialgleichungen und bei Differenzen- 
gleichungen). Es folgt ein umfangreiches Literatur- 
verzeichnis. Die äußere Form des Berichts ist recht 
ansprechend; verschiedene Schriftarten sorgen für 
Übersichtlichkeit, und Abkürzungen der häufigsten 
Begriffe gewährleisten eine knappe und prägnante 
Ausdrucksweise. 


Berlin R. Reıssıe 


Dr. K. Hirschfeld, Baustatik. Theorie und Bei- 
spiele. XVI + 8238. m. 1263 Abb. (1849 Einzeldar- 
stellungen) im Text und in 217 Zahlenbeispielen sowie 
38 Zahlentafeln m. 560 Fig. Berlin/Göttingen/Heidel- 
berg 1959. Springer-Verlag. Preis geb. 76,50 DM. 


Dieses Lehrbuch über Baustatik behandelt in den 
ersten fünf Kapiteln die Stabstatik. Es werden be- 
sprochen die Grundlagen der Statik, das Verfahren 
der Kinematik, Formänderung, Verfahren zur Auf- 
lösung linearer Gleichungen, Berechnung statisch 
unbestimmter Systeme nach verschiedenen heute 
gebräuchlichen Verfahren (Kraftgrößenverfahren, 
Formänderungsverfahren, Festpunktverfahren, Itera- 
tionsverfahren, Verfahren vonK ani, kombinierte Ver- 
fahren). Einen besonderen Hinweis verdient die ein- 
gehende Behandlung des Vierendeelträgers und eini- 
ger räumlicher Stabwerke. Ein weiteres Kapitel 
befaßt sich mit den Balken und der Kreisplatte auf 
elastischer Unterlage. Es finden sich dann noch eine 
Sammlung über Aufgaben und eine große Anzahl sehr 
wertvoller Hilfstafeln, die ein schnelles Auffinden von 
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Ergebnissen ermöglichen. Das Buch wird abge- 


schlossen mit einer Literaturschau und einem Schlag- 


ortverzeichnis. 
ee der sehr sorgfältigen Auswahl des Stoffes 
zeichnet sich das Werk vor allem durch die große An- 
zahl der durchgerechneten Beispiele aus, bei denen 
alle für die jeweilige Aufgabe erforderlichen Bechen- 
schritte auch tatsächlich ausgeführt sind. Der er- 
läuternde Text ist klar, aber nicht zu breit abgefaßt. 


Das Buch stellt zweifellos eine wertvolle Bereiche- 
rung des Schrifttums über Baustatik dar, insbesondere 
weil es die vielen unterschiedlichen Verfahren, die 
meist nur in verschiedenen Büchern zu finden sind, in 
einem Werk vereinigt. Es dürfte den Studierenden 
und dem in der Praxis tätigen Ingenieur eine wirk- 
liche Hilfe sein. 


Hannover W. ZERNA 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten). 


J. F. Holbrock, Laplace Transforms for Elec- 
tronie Engineers. XIII + 2598. m. 170 Abb. 
London 1959. Pergamon Press. Preis geb. 50 s net. 


V. Fock, The Theory of Space, Time, and 
Gravitation. XVIII + 411S. London 1959. Per- 
gamon Press. Preis geb. £ 5 net. 


6. Asser, Einführung in die Mathematische 
Logik, Teil I. IV + 1848. Leipzig 1959. B.G. 
Teubner — Verlagsgesellschaft. Preis geb. 11,25 DM. 


W. Olszak, Non-Homogeneity in Elastieity 
and Plastieity. VII + 5288. London 1959. Perga- 
mon Press. Preis geb. £ 5 net. 


P.A.P. Moran, The theory of storage. 1108. 
London 1959. Methuen & Co. Ltd. Preis geb. 
13 s 6 d.net. 


R. Sauer, Ingenieur-Mathematik, Erster 
Band. VIII + 304 S. m. 178 Abb. Berlin/Göttingen/ 
Heidelberg 1959. Springer-Verlag. Preis geb. 24,— DM. 


Zuschrift an den Herausgeber 


A note on my paper „A Geometrie Derivation of 
the Non-Linear Stress-Strain Relations for Isotropie 
Elastie solids*. 

In a previuos paper (ZAMM 39 (1959), p. 423—424) 
there was introduced an isotropie tensor Hijxı with the 
assumption „that the elasto-mechanical properties 
of an isotropie solid... .in E, are entirely characterized 
by the isotropie tensor Eijr1“*. Theinevitable question 
is: what is the nature of this — inuitively introduced 
r tensor ? We shall try there to give an interpretation 
of it. 

* 

Dince the difference Bijrı — Eijkı is of the dimen- 
sion of stress, it is evident that the dimension of Eijkı 
has to be the same. Hence Eijkı should be the fourth 
order tensor of the initial stress, due to the forces of 
cohesion acting between particles of the solid. When 
such forces are equal in all direetions, the field of 
such forces is to be represented by a tensor whose 
components are invariant with respect to the rotation 
of coordinate axes — and that is the tensor Eijkı. 

Moreover, when the deformation of a solid is 
affected the tensor Hijkı as a geometrie object defined 
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3. Horn/H. Wittieh, Gewöhnliche Differential- 
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C. €. Lin, Turbulent flows and heat transfer. 
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over particles of the solid deforms itself and the de- 


* 
formed tensor Eijrı describes the field of cohesion 
forces after the deformation, while the stress tensor 
tijkı presents the change of the cohesion stress due to 
external forces that caused the deformation. 
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 Biegungen und Schwingungen linearer Kontinua behandelt (gespannte Fäden mit diskreten Massenpunkten sowie Saiten | 

5 und. Stäbe mit ‚beliebiger Massenverteilung). Dabei werden die sogenannten Oszillationseigenschaften, die « die Einfach- 
Kr = der Eigenfrequenzen sowie die Anzahl und Lage der Knoten der Eigenschwingungen betreffen, eingehend diskutiert; 
Be sie hängen mit gewissen Vorzeicheneigenschaften der auftretenden Matrizen und Kerne eng zusammen, Die Theorie ist en 
Er u 2 diesem Buche erstmals vollständig UmEEoSNellE, dabei sind auch die zahlreichen einschlägigen Arbeiten der Are "= 
| übersichtlich zusammengefaßt. 


AL Für ‚den reinen Mathematiker und jeden, der mit solchen kn, ‚und ren zu tun har sinddie 
5 7 erörterten Probleme von Wert und Interesse. Hier ran: die klare und ee a den Bedürfnissen 

I 5 des Praktikers. ’ 

z In einem Anfangskapitel de die wesentlichen Sätze der Werisenähenzie le Daher genügen zum Verständnis 

er des Werkes Kenntnisse über Differential- und Integralrechnung. Lediglich die letzten nn setzen die Kenntnis linearer 
Integralgleichungen voraus. : 


” AUS DEM INHALT 
Kapitel I Grundlegendes über Matrizen und quadratische Formen 


I KapitelII _Oszillationsmatrizen 
 _ KapitelIII- Kleine Schwingungen mechanischer Systeme mit n Freiheitsgraden 


5 KapitelIV Kleine Schwingungen mechanischer Systeme mit unendlich vielen Freiheitsgraden ! 

7 , 2 4 DER En 
I KapitelV _Zeichenfeste Matrizen a 
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Bestellungen durch eine Buchhandlung erbeten 


ATZE A-D.- EM TE». YVVE RT 


In dem Werk sind alle Te ya 
DDR im November 1958 in Jena ak ing Ai 3 
enthalten. Durch die Beteiligung auch Br 
erhält das Buch internationale Bedeutung. Hohandeh me A 
forschung, wie ‚‚Über eine Verallgemeinerung der 
und Ausführungen über „Optik, Physik und physikalische € 
italienischer Sprache. Neben Problemen der Halbleiter werden F ze kt 
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Schwingungen aller Wellenlängen des Erdkörpers besprochen. 


